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This paper seeks to demonstrate that many of the existing mass-measurement variables proposed 
for hadron coUiders (mr, ?Tieff, wit2, missing pr, hr, v'^min, etc.) are far more closely related 
to each other than is widely appreciated, and indeed can all be viewed as a common mass bound 
specialized for a variety of purposes. A consequence of this is that one may understand better the 
strengths and weaknesses of each variable, and the circumstances in which each can be used to best 
effect. In order to achieve this, we find it necessary first to revisit the seemingly empty and infertile 
wilderness populated by the subscript "T" (as in "pr") in order to remind ourselves what this 
process of transversification actually means. We note that, far from being simple, transversification 
can mean quite different things to different people. Those readers who manage to battle through 
the barrage of transverse notation distinguishing "T" from "V" or from "o", and "early projection" 
from "late projection" , will find their efforts rewarded towards the end of the paper with (i) a better 
understanding of how collider mass variables fit together, (ii) an appreciation of how these variables 
could be generalized to search for things more complicated than supersymmetry, (iii) will depart 
with an aversion to thoughtless or nai've use of the so-called "transverse" methods of any of the 
popular computer Lorentz-vector libraries, and (iv) will take care in their subsequent papers to be 
explicit about which of the 61 identified variants of the "transverse mass" they are employing. 
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I. INTRODUCTION 

Almost every analysis of data from hadron colliders 
uses at some point a variable which represents a "pro- 
jection" of an energy or momentum into the plane trans- 
verse to the beams. The typical reason for performing 
these projections is that one does not wish the analysis 
to be sensitive to the unknown momentum - along the 
direction of the beams - of the quarks or gluons which col- 
lide in the 'hard' interaction. Given the widespread use 
of such variables it is perhaps surprising that many col- 
lider physicists are probably unaware that there exist at 
least two commonly-used ways of projecting of a Lorentz 
energy-momentum vector into the transverse plane, and 
that these two different methods have very different prop- 
erties when the mass is non-zero (see Section III below). 
Furthermore, as explained later in Section [V] for each of 
those transverse projections, there are at least two in- 
equivalent ways that transverse vectors can be "added 
together", each of which has benefits and weaknesses. A 
careful definition of what we mean by a transverse pro- 
jection forms the first part of this paper. 



The later part of the paper (Sections VI -XI I deals 
with mass-scale (or energy-scale) variables, a variety of 
which have been proposed in the run-up to the LHC data- 
taking^ . Though some of these variables have been con- 
structed from careful consideration of the Lorentz sym- 
metries of space-time, others have been created in a some- 
what ad-hoc process, after simulations demonstrate that 
they provide good signal-to-background discrimination, 
or that they are highly correlated with the mass of some 
particle or particles. The main aim of this part of the 
paper is to demonstrate that many of these seemingly ad- 
hoc definitions are in fact not only well-motivated from 
the kinematical perspective, but also that the associated 
variables are more closely related than one might have 
thought. 

Figure[T]illustrates some of the variables that are found 
to be connected in ways that are not widely appreciated. 
One might argue that we add little to the sum total of 
human knowledge by merely showing the relationships 
between existing variables which are already known to 
work well in particular roles. However, careful study of 
their similarities and differences not only gives insights 
into why (and under what circumstances) these choices 
are appropriate, it also fits them into a common frame- 
work - from which it is straightforward to make general- 
izations to more complex decay topologies. 

The paper is organized as follows; first we carefully 
define our notation for Lorentz 1+3 vectors and their 
transverse projections in Section |llj Then in Section |III| 
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FIG. 1: The stretched, webbed limbs of the Glaucomys volans 
have been adapted by generations of natural selection to pro- 
vide an ideal visual illustration of the various different, yet 
related, transverse mass variables (and incidentally provide 
an appropriate aerodynamic shape for gliding flight). 

Photograph © Joe McDonald. 



we describe the two common but inequivalent transverse 
projections, which we shall denote by subscripts T or V. 
We also introduce the special case of a "massless" trans- 
verse projection, denoted by o. In Sect ion [TV] we compare 
the results from the three different types of projections: 
T , V and o . In Section |V] we highlight the differences 
between projecting into the transverse plane before or af- 
ter forming composite objects. Section |VI| describes the 
general event topology targeted by new physics searches 
in channels with missing momentum. 



^ For a recent review see [l]. 



All of those ingredients are put to work in Sections VII 
to |Xj which contain the main results of this paper. 
In Section |VII| we introduce the general class of mass- 
constraining variables which can be usefully applied for 
studying events containing invisible particles. The set 
of possible transverse mass variables is extended in Sec- 
tion [Vlllj where we consider additionally projected one- 
dimensional objects. Some mathematical properties of 
these mass-constraining variables are discussed in Sec- 
tion |IX| Some of the variables have previously appeared 
elsewhere in the literature and we clarify the correspond- 
ing connections in Section |Xj In Section pQ| we illustrate 
the use of these variables with two simple examples: an 
s-channel resonant production process, for which we take 
inclusive Higgs boson production pp h W^W~ — >■ 
£^£^ -t- frpi ^nd a pair-production process represented 
by top quark production pp ^ tt ^ hM^i~ -\- $rp. Sec- 
tion [XlT] contains a short summary and conclusions. 

Appendix |A] contains a short guide to the currently ex- 
isting computer libraries and codes which can be used for 
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computing some of the variables described in the main 
body of the text. Appendix [B] provides derivations of 
extremal mass-bound results and other general mathe- 
matical proofs which are used elsewhere in the paper. 



II. NOTATION AND CONVENTIONS 

A. Labelling momenta and their components 

In general, capital letters {P, Q, AI, E, etc.) will refer 
to genuine 1-1-3 dimensional vectors, while lowercase let- 
ters {p, q, m, e, etc.) will refer to "less than 1-1-3" dimen- 
sional constructs. Lower indices . . . label individual 
final state particles, while lower indices a,b, . . . are used 
for parent particles and the corresponding collections of 



final state particles defined below in Sec. |VI| We also use 
upper indices fj,,v, ... to label the components of 1-1-3 vec- 
tors, and upper indices a, /?, . . . to label the components 
of the projected 1-1-2 dimensional transverse "vectors" of 



the types defined in Section III The 1-1-3 metric g^jy is 
diag(l, — 1, — 1, — 1) and the 1+2 dimensional metric g^^ 
is diag(l, —1,-1). Thus the 1+3 energy-momentum vec- 
tor for some particle is written P'^ — {e, and the 
corresponding mass denoted by a capital M: 



M' 



P\ (1) 

As illustrated in Fig. [2] any 3-dimensional vector P 
can be trivially decomposed into a transverse and a lon- 
gitudinal component: 



P=ipT,Pz)- 



(2) 



The transverse momentum px — {px,Py) of the parti- 
cle is, of course, 2-dimensional, so it has a lowercase 
"p" . Similarly, the longitudinal momentum pz is 1- 
dimensional, and is also lowercase. By contrast, the en- 
ergy E measured in the detector is a component of a "1+3 
dimensional thing" , since it is given in terms of the 1+3 
dimensional mass AI and the 3-dimensional momentum 
P: 



E = \ AP 



)2 _ 



+ Pt+Pz 



(3) 



When it comes to projecting geometric 3-vectors like 
P, the decomposition shown in eq. ^ and Fig. [2] is un- 
ambiguous. One has no other choice - the very definition 
of the transverse plane requires one simply to dispose of 
the z-component to arrive at pr = {Px,Py)- AH the trans- 
verse projections considered in this paper (and any others 
that one might invent) must share this property, or else 
they cannot justify being so named. 

However, "projecting" the time-like component E is 
not, in itself, a well defined operation. What does it 
mean? There is not a single correct answer, but rather a 
number of different answers, each with different proper- 
ties and motivations. How one should (and even whether 
one should) project time-like components of 1+3 Lorentz 
vectors is dependent on what one is trying to achieve. 
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FIG. 2; The standard geometry of a collider experiment. The 
z axis (in blue) is oriented along the beam, while the x and y 
axes (in red) define the transverse plane. Any 3-dimensional 
vector P can be uniquely decomposed into a longitudinal com- 
ponent pz and a transverse component pr. 



B. Labelling transverse projections 

In the particle physics literature, one can find evidence 
of at least three different types of "transverse projection" 
being applied to (l+3)-Lorentz vectors - although this 
diversity is not obvious at first glance, as the majority of 
papers do not explicitly state which projection they are 
using. ^ Even those papers which define the projection 
explicitly, usually neither comment on why the particular 
choice was made, nor comment on what would happen 
were another projection to have been used. 

One of the main objectives of this paper is to place 
these three main types of projection side-by-side so that 
their differences, the things that they share, and their 
respective uses can be directly compared. Before we de- 
scribe them in more detail, we make some remarks about 
notation. 

In the literature, all of the types of projection are la- 
belled by the same symbol: the letter "T" . Since in this 
document we need to clearly distinguish the three types 
of projection, it is necessary for us to create our own no- 
tation for each - and we use the three symbols "T" , "V" 
and "o" for that purpose. 

We will continue to use the letter "T" to indicate 
"generic" transverse quantities, i.e. quantities which are 
either common to all projections (e.g. the transverse 
momentum (2)-vector already commented upon, the 
missing tranverse momentum vector fj,, or the transverse 
upstream visible momentum vector ut defined below in 



Section VI A I or for quantities which for historical rea- 
sons carry a transverse subscript, but which may not be 
tied to one type of projection to the exclusion of others 
(e.g. Ht). 

Note that certain quantities, such as the so-called 



^ This may be because all forms turn out to be equivalent for 
massless particles. 
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"transverse energy" and "transverse mass", are differ- 
ent in each of the projections. For this reason the sym- 
bol "ct" is effectively meaningless, and should appear 
nowhere in this document (outside this sentence) un- 
like ct, Cv and Co (which are all different and all well- 
defined). Similarly, rriT is also ambiguous, and should 
be specified as being m-j, my or ruo- In contrast, "pr" 
is perfectly legitimate, and indeed (as we have already 
noted) is equivalent to pV , Pv and po : 



III. TRANSVERSE PROJECTIONS 



(4) 



In this Section we describe the three different types of 
projections "T", "V" and "o". While reading this and 
the following sections, the reader may find it helpful to 
refer to Table |l] for notational reference, and also to see 
how the results for each projection compare to those of 
the others. 



A. The mass-preserving "T" projection 

The first approach we will describe, which will be de- 
noted by a "T" subscript, is the most common in the 
mass measurement literature. For example it is found in 
the early literature on the transverse mass when it was 
used to measure the W mass [2H5] and in the generaliza- 
tion of the transverse mass to pair production, namely 
Mt2 (the stransverse mass) [71-1^ as well as in literature 
relating to Mct ;30-3S] and in reviews of the field [1 . 

In the T projection one defines the 1-1-2 dimensional 
transverse energy^ ct and transverse momentum pV in 
terms of the l-t-3 dimensional mass M and 1-1-3 dimen- 
sional components according to 



my = M. 



(5) 
(6) 
(7) 



In this case, the components of the 1+2 dimensional 
quantity 



Pt = (eT,PT) 



(8) 



^ Note that in equation (jsj, it is the middle expression + p^, 

that we use to justify our "calling" the LHS a (transverse) "en- 
ergy" - since it is square root of a "mass squared plus a trans- 
verse momentum squared" . Someone who saw the right hand 
expression first, E'^ — p| , could argue differently, and might 
reasonably expect us to call the whole quantity a "longitudinal 
mass" - since it is a square root of an "energy squared minus a 
longitudinal momentum squared" . All this really goes to show 
is that the "name" of the quantity is to some extent a matter of 
convention rather than physics. 



satisfy the mass shell condition 



(9) 



with the 1-f 3 dimensional mass M. 

The equivalence class for this projection function - the 
set of 1-1-3 vectors which map to the same 1+2 projected 
vector under T - consists of the set of 1-1-3 vectors with 
the same pV and M: 



A-P+pI,Pt 



(10) 

The fact that all members of the equivalence class share 
the same mass is what motivates us to call this the "mass 
preserving" T projection. 

Given its dominant use in the literature, it is something 
of a surprise that the nomenclature of the T projection 
is not adopted in the commonly used high-energy physics 
computer libraries such as CLHEP PH] or ROOT [37] which 
instead implement the alternative V projection intro- 



duced below in Section [III B[ The T projection is, how- 
ever, used in the "Oxbridge stransverse mass library" [38] 
and the U.C. Davis Mt2 library [3S]. See Appendix [AI 



and Table VII in it for a summary of library conventions. 



B. The speed-preserving "V" projection 

Alternatively one can follow the method of the CLHEP 
[36] and ROOT [37] libraries and "project" the energy on 
the transverse plane, using the same angle 9 as for the 
momentum vector. As alteady seen in Fig. [2] the mag- 
nitude Pt of the transverse momentum pT is related to 
the magnitude P of the 3-dimcnsional momentum P by 



Pt = P smt 



with 



tan 9 ; 



Pt 

Pz ' 



(11) 



(12) 



Thus by analogy with ( 11 ) one can define the transverse 



energy in terms of its 1+3 dimensional counterpart E as 

ev = Ssin6'. (13) 

Then for any individual 1+3 momentum vector we have 
the V version of the "transverse" components 



e\j — E sin 9 = 



Pt 



Pv 
m 



V 



Pt, 

M sin e 



Vpt+pI 

Pt 

Vpt + pI 



E, 



M. 



(14) 

(15) 
(16) 



We can take the angle 9 to be defined in (0, tt), so that 
ev and TOv are always nonnegative. 



6 



Quantity 


Transverse projection method 


Mass-preserving 'T' 


Speed-preserving 'V' 


Massless 'o' 


Original (4)-momentum 
(H-3)-mass invariant 
Transverse momentum 


M = ^E''-p}-pi 

PT = {Px,Py) 


(H-2)-vectors 


Pt = (eT,PT) 


Pv = (ev,Pv) 


Po = (eo,Po) 


Transverse momentum 
under tlie projection 


PT = PT 


pv = PT 


Po = PT 


Transverse energy 
under tlie projection 


ex = y/M^+p} 


ev = E\sine\ = \pt\/V 


Co = \pt\ 


Transverse mass 
under tlie projection 


2 2-2 

mj- = ey -px 


2 2 ^2 

mv = ev - Pv 


2 2 -2 r\ 


Relationship between 
transverse quantity and its 
(1+3) analogue 


my = M 


mv = M \sm9\ 


ruo — 






Wo = 1 


Equivalence classes under 
(1 + 3) ^ (1 + 2) 


All P*^ with the same 
Px, Py and M 


All P^ with the same 
Px, Py and V 


All with the same 
Px and Py 



TABLE I: A comparison of the three transversification methods introduced in Section [Till 



In this V method of projection we can also introduce 
1+2 "vectors" which now have components 



K = (ev,Pv) 



(17) 



The V projected components obey a different mass sheU 
relation than the T projected components in ([9|: 

el-pl=ml<M\ (18) 

with the 1+2 dimensional V projected mass my. 

Just as an aside, one could also define the "longitudi- 



nal" components in complete analogy to (14)-(16l 



Vz 



= £;|cos6i| = 

= Pz, 

= Af I COS 6*1 = 



\Pz\ 

Vpt+pI 



(19) 
(20) 
(21) 



although in what follows we shall not be making any use 
of those. The connection between the 1+3 dimensional 
quantities and the V 1+2 dimensional components is 



2 , 2 

ew + e , 



(22) 
(23) 

For massive vectors'' the equivalence classes of the V 
projection are different from those of the T projection. 



The mass-shell relation (18) implies that all the 1+3 vec- 



tors which map to the same 1+2 vector under the V pro- 
jection share the same value of my = M sinO and thus 
generally do not preserve the usual invariant mass M, 
since my ^ M for any ^ f . 

A more physical picture of the equivalence class of vec- 
tors for the V projection can be found by considering the 
3-speed of the particle 



V = 



P 
E' 



(24) 



After the V projection, the corresponding 2-speed is given 
by 



Vy 



Py Pt PsinO P 
Cy Cy E sin 9 E 



(25) 



Eqs. (24 1 and (25) reveal that the V projection is "speed 



preservmg , i.e. 



Vy = V, 



(26) 



which justifies our choice of subscript notation for this 
kind of transverse projection. The equivalence class for 
the V projection therefore consists of all 1+3 vectors with 
the same px and speed V: 



Vpf 



Pz 



V 



,PT:Pz 



PT 

V 



,PT 



(27) 



* See section 



IV B 



for comments concerning the massless case. 



Note that members belonging to the same equivalence 
class under the V projection (27) have the same speed, 



7 



but different masses, while members of the same equiv- 



alence class under the T projection ( 10 1 have the same 
mass, but different speeds. 



The massless 



projection 



The massless "o" projection defines components 



eo 

Po 



\Pt\ 
Pt 



(28) 
(29) 



and thereby defines a massless 1+2 vector of the form 
Po^{\ptIpt). (30) 

The main feature of this projection is that the 1+2 vector 
always has a null invariant 



gapPoPo =ml=0. 



(31) 



It should be noted that p" and p" have three degrees of 
freedom {{bt ,Px,Py} and {e\/,px,Py}, correspondingly). 
Therefore their equivalence classes are one-dimensional, 
and can be parameterized by the coordinate Pz, as in- 
dicated in (10) and (27 1. In contrast, our 'o' projected 
vector p" has only two degrees of freedom, p^ and Py — 
the time- like component being fully specified from p^ and 
Py through Co = \pt\- The equivalence class of any p" 
vector is therefore also a 4 — 2 = 2-dimensional object, 
parameterized by, say, Pz and E: 



{E,PT,Pz) 1^ {\pt\ ,Pt) ■ 



(32) 



IV. 



COMPARISON OF THE DIFFERENT 
TRANSVERSE PROJECTIONS 



The three different projections discussed in Section [TTT| 
are pictorially represented in Fig. [3j For a given fixed 
value of Pt, the white region in the figure depicts all pos- 
sible allowed values of the energy E and the longitudinal 
momentum pz ■ (The yellow-shaded region E^ < p^ + p^ 
is forbidden because it corresponds to a tachyonic par- 
ticle with < 0, travelling with superluminal speed.) 
In this figure, we consider the plane of energy squared 
versus momentum squared, and in order to retain the in- 
formation about the sign of the longitudinal momentum 
component, we plot sigTi{pz)p1, so that the mapping from 
the {E,pz)-p\'Ane to the (£'^, sign (p^) p^) -plane is one-to- 
one. 

Each of the three transverse projections maps a point 
with some given^ values of E and Pz onto the pz — 
axis as shown. In the case of T, the projection is along a 



For definiteness, in Fig. [s] we have chosen an illustration point 
with Pz < 0. 




FIG. 3: A pictorial representation of the three transverse 
projections discussed in Section [HI] The colored arrows rep- 
resent the mappings under the T, V and o projections. The 
blue and green dotted lines represent the equivalence classes 
of the projected points under the T and V projections, re- 
spectively. 



line of constant mass M and results in transverse energy 
squared = + Pt- In the (i?^, sign(p^) p^)-plane, 
lines of constant M are straight lines, which explains our 
choice of quadratic power scale on the axes. Fig. [3] il- 
lustrates that the equivalence class of vectors under the 
T projection is one-dimensional: it is represented by the 
two blue dotted straight lines, which can be simply pa- 
rameterized by the value of Pz ■ 

The V projection, on the other hand, projects along 
a line of constant speed V, as indicated in Fig. [3] In 
the {E'^ ,sign{pz) pl)-plan.e, lines of constant V are also 
straight lines, albeit with a different slope. The result- 
ing value of the transverse energy is Cv = pt/V. The 
corresponding equivalence class of vectors is given by the 
two green dotted lines, and can also be parameterized in 
terms of a single parameter, say Pz- 

Finally, the massless "o" projection maps any allowed 
point in the {E,pz)-plane to the massless 1+2 vector with 
transverse energy Co = Pt- The equivalence class of vec- 
tors in this case is two-dimensional, and is represented 
by the whole white shaded region in Fig. [3j 

All of the previous discussion can be recast in the lan- 
guage of the {A'P,sign{pz)pl)-p\ane, as shown in Fig. [Z] 
In this case, the whole > half-plane is allowed, and 
the T projection projects horizontally onto the Pz — 
axis, following the blue arrow. The V projection is also 
done along a straight line, following the green arrow. As 
before, the equivalence classes for the T and V operations 
are straight lines, while the equivalence class for the "o" 
case is given by the whole > half-plane. 
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FIG. 4: The same as Fig. |3] but plotted in the 
(M^,sign(p2)p^)-plane. 



A. A hierarchy among projections 

As illustrated in Fig.|4j the definition of each projection 
imposes a hierarchy among the projected masses of the 
form: 



M — Tot > TOv > Too = 0. 



(33) 



We draw attention to this hierarchy here as it will have 
very close analogues in the more complicated experimen- 
tal mass bounds derived from each type of projection in 
the later sections. 

Given the mass hierarchy ( [33| , eqs. (|9]), (18) and (28) 
imply that a similar hierarchy exists for the projected 
energies: 

> ct > ev > Bo = pt, (34) 
which is illustrated in Fig. [3] 

B. Equivalence in the massless Hmit 

We note that in the special (but common) case in which 
the original four-vector is massless (M = 0) all projec- 
tions are equivalent since 



lim ct 



lim 

M-i-O 



\pt\ 



(35) 



and thus 



lim - lim K =p: = {\ptIpt) ■ (36) 

Clearly the projections are not equivalent for massive 
particles, nor for collections of massless particles (unless 
they be collinear) since collections of massless particles 
can have large total invariant mass - the equivalence ex- 
tends only to application to individual massless particles. 



In practice, the statement above may also be taken as 
saying that all the projections are equivalent in the high- 
energy limit - i.e. the limit in which the momentum of 
a particle is much greater than its mass - again only at 
the level of individual high-energy particles. 

Since all the projections are equivalent in the above 
limits, and since most individual reconstructed particles 
in high-energy physics experiments satisfy one of those 
limits due to the small masses of the leptons and light 
quarks, one might wonder what all the fuss is about. 
However, the importance of the distinctions will be seen 
to arise and become very large when we consider com- 
posite particles, i.e. collections of massless "daughter" 
particles^ . Composite particles are expected to have non- 
negligible masses, even when they consist of sums of (ap- 
proximately) massless particles. As we already learned 
from the simple example considered in Figs. |3]and|4) not 
only do these composite particles generate very differ- 
ent projected l-t-2 vectors, but the classes of equivalent 
four-vectors associated with those projections are very 
different as well. 



V. SUMMING AND PROJECTING: EARLY 
VERSUS LATE PROJECTIONS 

In forming transverse kinematic variables for compos- 
ite particles, one needs to perform two separate opera- 
tions: summation of the momentum vectors of the daugh- 
ter particles, and projecting into the transverse plane. 
The order of these operations does not matter for the 
two space-like vector components: 



i 
i 




(37) 
(38) 
(39) 



where we use an index i to label the momenta of the 
individual daughter particles and the sums run over all 
such daughter particles^. 



® The need for considering composite particles arises when deahng 
with short-lived heavy resonances, which decay promptly to a 
certain collection of daughter particles, which in turn are seen in 
the detector. The energy and momentum of the parent resonance 
are correspondingly obtained by summing the measured energies 
and momenta of the daughter particles. 

Recall our convention that lowercase letters refer to 1+2 dimen- 
sional quantities and capital letters refer to 1-1-3 dimensional 
quantities. Thus in the left-hand-sides of eqs. (j37]l-|39]l we are 
adding 2-dimensional transverse vectors, while in the right-hand- 
sides we are first adding the corresponding 3-vectors, then pro- 
jecting their sum onto the transverse plane. 
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However, projecting before or after the sum can make 
a very significant difference to the value of the time-hke 
(ey, Cv or Co) component of the final 1+2 vector - and 
therefore the operations of projecting and summing do 
not generally commute: 



E 




(40) 
(41) 
(42) 



One can see clearly how the order makes a difference 
if one considers an extreme case consisting of a pair of 
massless daughter particles travelling in opposite direc- 
tions along the beam pipe, i.e. with 1+3 momenta 



Pf = (£^,0,0,+^;), 
= {E,0,0,-E). 



(43) 
(44) 



If one were to project these 1+3 momenta into the trans- 
verse plane before summing (a combined operation here- 
after called early projection) , one would find that the re- 
sulting 1+2 dimensional vector 



T.P'^ = P1T+P2T 



= iE,0,0,E)^ + {E,0,0, 
= (0,0,0) + (0,0,0) 
= (0,0,0) 



-Eh 



(45a) 

(45b) 
(45c) 
(45d) 



is null. A null sum would also be obtained if we had used 
the V or o projections.® However if one were first to sum 
the Lorentz 1+3 vectors P/" and then later project into 
the transverse plane (hereafter denoted late projection) 
one would find that 



{Pi: 



(46a) 



= HE, 0,0, E) + {E, 0,0, ^E))^{46h) 
= {2E,0,0,0)^ (46c) 
= {2E,0,0), (46d) 



which is clearly not the same as was found in ( 45d ) . This 



extreme case shows that while projecting early has the 
effect of reducing dependence on longitudinal momenta, 
projecting late means that the resultant projected com- 
posite retains much more sensitivity to the original rela- 
tive momenta along the beam directions. 



* In fact for this example we have chosen massless vectors for which 
the 'T', 'V', and 'o' projections are identical. 



This concludes this section, whose main purpose was 
simply to highlight the difference between the "early" 
and the "late" transverse projection. It also underscores 
the need to develop the proper notation to distinguish be- 
tween these two types of transverse projections, which we 
shall do below in Section fVI CI The differences between 
the two projections will be further illustrated with the 
physics examples considered in the later sections. One 
may reasonably wonder which one of the two projections 
is more appropriate and should be used. In principle, 
the answer to this question will depend on the analysis 
being performed. If one is initially building a compos- 
ite particle from two leptons, e.g. from a Z-boson decay 
Z — > e^e^, then the relative longitudinal momentum of 
the positron and the electron is probably a safe quantity 
to retain full sensitivity to in one's calculations. However, 
in cases where jets at large rapidity \r]\ are concerned, the 
probability of QCD radiation grows rapidly as one gets 
closer and closer to the beam direction. One will often 
prefer not to have the high-energy end of the composite- 
particle spectrum dominated by combinations of low \pT \ , 
high-energy forward-going jets with other low \pt\, high- 
energy backward-going jets, so in this latter case, early 
projection would probably be appropriate. Nevertheless, 
giving a universal prescription for selecting the "correct" 
transverse projection for collections of particles is beyond 
the scope of this paper. The best method will depend 
on non- kinematic factors, such as the size of any back- 
grounds, the detector resolution, and other factors that 
will vary from case to case. 



VI. INTERPRETING EVENTS 
A. Characterizing an event 

Analysis of an event is a game. The aim of the game 
is to interpret the available information within a par- 
ticular framework or hypothesis. In this paper we wish 
to employ a very general framework that will be useful 
for searches and mass measurements at hadron colliders 
{pp, pp or even pp for that matter). Specializations of 
this framework will then be useful in a wide variety of 
different contexts. The general layout of an event is rep- 
resented in Fig. [5] The figure comprises: two incoming 
objects, denoted by the proton lines on the left hand side; 
an interaction, represented by an oval 'blob'; and some fi- 
nal state objects, contained within the rectangles on the 
right hand side. Since it is the final state objects that 
provide the kinematic information about the event, we 
now take some time to explain rather carefully what we 
mean by them. 

We define final state objects of two types. A visible fi- 
nal state object is one that leaves a signal in the detector 
that betrays its presence. Those signals may then be re- 
constructed and interpreted as an individual particle - for 
example as photon, electron or muon - or the signals may 
be indicative of a composite object, such as a QCD or tau 
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visible momentum 
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Visible daughters Vi 


Parent Pa 
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Invisible daughters I2 
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Visible daughters Vat 




Invisible daughters In 



FIG. 5; The event topology for new physics searches and 
measurements used in this paper. 



jet. The "visible object" category is deliberately allowed 
to be sufficiently broad as to permit the inclusion of very 
heavy, visibly decaying, composite objects such as Z, W 
or top quarks. The classification of a final state object 
as 'visible' here implies not only that a signal consistent 
with the presence of some particle has been observed, 
but also that the full Lorentz energy-momentum vector 
of that particle can be reconstructed from the observed 
signal (to within some experimental precision). For most 
heavy visible objects (jets, Z, H bosons, ...), the 
four-momentum of the visible object must be calculated 
from the vector sum of its constituents. 

By contrast, an invisible final state object is one that 
leaves no direct signal, but the existence of which is de- 
manded by the interpretation of the event being imposed 
by the analyst. The numbers, types and masses of any in- 
visible final state particles form part of the interpretation 
of the event. The 3-monicntuni vectors of all invisible 
particles are a priori unknown, and are constrained only 
by conservation of the total momentum of the event in 
the plane transverse to the beam. The general framework 
can accommodate a final state hypothesis in which invis- 
ible particles have particular known (or rather assumed) 
masses, but it can also be applied when some or indeed 
all of those invisible particles have unknown masses. 

As illustrated in Fig. [5] the next step in interpreting 
the event is to partition the combined set of all final state 
objects (visible and invisible) into subsets, which are rep- 
resented by rectangles in the figure. Each final state 
object must be found in one and only one such subset. 
There is one subset per parent plus one further subset, 
the latter being labelled "upstream visible momentum" 
in the figure. 

In our interpretation a parent is any short-lived ob- 
ject that is believed to have decayed to produce the vis- 
ible and invisible final state objects in its associated set. 
The term 'parent' is usually associated with a short-lived 



heavy state, most often a reasonably narrow resonance 
(whether produced directly in the "hard scatter" or from 
decays of even heavier objects). 

The general framework presented permits a variety of 
different interpretations for any given event. For any 
particular interpretation there is a corresponding par- 
titioning of the final state into subsets. For example 
an event which contains evidence of an electron and a 
positron, and which is hypothesized to also contain a 
neutrino and an anti-neutrino, could be partitioned into 
parent/daughter combinations: — )■ {e"*", Ve} and 
{e^, Ve] for one analysis; however another anal- 
ysis might find it more appropriate to partition those 
objects according to the interpretation — > {e+, e~} 
and Z° {ve, Ve\- 

The subset corresponding to any parent may contain 
any number (including zero) of visible particles and any 
number (including zero) of invisible particles — though 
it is not meaningful to have a totally empty set of daugh- 
ters. The framework is very general, in that the number 
of parents can be arbitrary, and the nature, mass and 
decay mode of any parent need not be related to those 
of any other. There is therefore a great deal of freedom 
in performing the partition into subsets. We shall later 
be constraining the masses of the parents so the subsets 
should be chosen to correspond to the descendants of the 
parents whose invariant masses we are interested in. 

Figure[5]also shows the one special (non-'parenf) sub- 
set into which visible final state particles may be allo- 
cated. That set is labelled ^^upstream visible momentum 
(UVM)", and is designed to be a 'catch-all' that will ac- 
commodate any visible particle not allocated to any of 
the parent sets. This is a special set in the following 
senses: firstly it is permitted to contain (by assertion) 
only visible objects; and secondly, and crucially, final 
state objects allocated to this set are not used directly to 
constrain the mass of any parent. Objects in this set are 
only used to keep track of overall energy-momentum con- 
servation. We do not specify the elements found in this 
UVM set, but in practical applications it almost always 
contains some contribution from "soft" particles that are 
unallocated to any parent. Such soft components usually 
include calorimeter energy found outside of jets, and low 
energy jets from multiple parton interactions, and per- 
haps from initial state radiation (ISR). As well as these 
soft components, one must include any other visible ob- 
jects not associated with any parent. The UVM set will 
often contain more than just 'soft' activity — since any 
type of visible particle can end up therein — possibly 
including decay products (of heavy progenitor particles) 
that the analyst chose not to allocate to any parent. In 
practice, every hadron collider event has some amount 
of UVM. Furthermore, as discussed in [101 US EH 121- 
[2ni 131] J the presence of a significant amount of UVM can 
in fact be beneficial in mass reconstruction studies. 

Apart from the reconstructed physics objects, another 
important experimental quantity is the missing trans- 
verse momentum in the event. This quantity is the exper- 
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imental collaboration's best estimate of the amount (and 
direction) of momentum in any particular event that has 
been carried away in the plane transverse to the beam 
by invisible particles. It is an important quantity insofar 
as we will wish to apply the constraint that the missing 
momentum in an event is entirely due to the invisible 
final state objects. 



B. Notation used to characterize events 

We require considerable amount of notation to describe 
events and the hypotheses and interpretations that we 
layer on top of them. We have summarized the nota- 
tion we have adopted in Table [ll| — and we recommend 
that readers immediately compare the first section of that 
table with any of the three small concrete examples pro- 
vided in Figures |6) [7] and [8] in order to follow later sec- 
tions. For the simplest pieces of notation, Table [IT| serves 
as the primary definition. Notation that requires more 
explanation will be described in more detail either below 
or at first point of use. 

The N parents are labelled Pq, (a = 1, 2, . . . , A^). The 
set of observed visible (hypothesized invisible) daughters 
associated with Pq is labelled Va (la)- Since no visible 
or invisible particle has more than one parent, we have 
^aCl^b = and Pi If, = when a ^ b, and so the num- 
ber of visible (invisible) particles may either be written 
as the sum of the number of visible (invisible) daughters 

of each parent Nv = Y.a=i {Nx = Y.a=i I^qI) or as 

the number of elements from the set of all visible (invis- 
ible) daughters TVy = |V| (iVi = where V = Uf=i Va 

As seen in Table [TTl in our conventions the letter "P" 
( "p" ) will be used to denote measured momenta, and the 
letter "Q" ("q") will be used for the momenta of any 
invisible or hypothesized particles. Correspondingly, the 
individual 4-momenta P/*, (i g V), of the visible daugh- 
ters are measured and known, while the individual 4- 
momenta , {i G I), of the invisible daughters are not 
measured and remain unknown. We denote the masses 
of the visible final state particles by Mi and those of the 
hypothesized invisible final state particles by Mi. Simi- 
larly, we will find it convenient to denote the 3-speeds of 
the visible final state particles as Vi and the 3-speeds of 
the hypothesized invisible final state particles by Vi. In 
some places we will need to refer to sets of these masses 
or speeds, and so we define: (i) the set consisting of the 
hypothesized masses of all invisible particles: 

tl = [M\iel}, (47) 

(ii) the set containing only the hypothesized masses of 
the invisible particles assocated with parent P^: 

fia = {M\iela} , (48) 



(iii) the set consisting of the hypothesized 3-speeds of all 
invisible particles: 

ij={v,\tel}, (49) 

and (iv) the set containing only the hypothesized 3- 
speeds of the invisible particles assocated with parent 
Pa: 

Va = {V,\tela}. (50) 

We denote the missing transverse momentum two- 
vector by the symbol^ and its magnitude thus p^. 

Note that some authors use variants of the symbol "-^7-" 
to denote the missing transverse momentum, but the 
distinction is necessary in this paper as we shall (as oth- 
ers should) make important distinctions between energy 
and momentum. 

We wish to apply the constraint that the missing mo- 
mentum in an event is entirely due to the Nx invisible 
particles with momenta , rather than to jet mismea- 
surement, for example. In other words, we use the rela- 
tionships expressed in: 

Nx Nv 

^ QiT = #r = -^PiT- (51) 

1=1 1=1 

in which the first equality represents our desire to con- 
strain the momenta of the invisible particles (and only 
those particles) using -/^j,, while the second equality re- 
minds us of our assumptions of how fj, is constructed as 
an experimentally measurable quantity. These relation- 
ships also remind us that we have assumed (i) that there 
are no sources of invisible momentum other than those 
coming from the parent decays, and (ii) that we have de- 
fined the "Upstream visible momentum" to contain all 
visible momentum deposits which did not originate from 
the decay of any parent. 

When considering the decay of a single parent P^ 

Fa^VaUla. (52) 



^ Note that due to its status as an experimentally measurable 
quantity, for the missing transverse momentum ^j, we use the 
letter "p" as opposed to "q", even though at high values ^j, is 
interpreted as the total transverse momentum of invisible parti- 
cles. 

" By right, since its meaning is derived from conservation of mo- 
mentum in the transverse plane, the missing transverse momen- 
tum ought universally to be known as fij,. Alas, much of the 
hadron-coUider literature, especially that from the experimental 
collaborations, calls the missing transverse momentum the "miss- 
ing energy" or "missing transverse energy" and denotes its mag- 
nitude "-^ji" and its two vector by some variant of 

"Et"- This 

is perhaps a result of history (a hang over from e"'"e~ or LEP 
terminology where the collision of point-particles from mono- 
energetic beams meant that one really could talk about missing 
energy) and the fact that is often reconstructed, at least in 
part, from calorimetric energy deposits under the assumption 
they were produced by massless physics objects. 
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Symbol 


Meaning 


See also 




1 A 1 


Cardmal number (number ot elements) ot any iimnte set A. 






Pa 


a**" parent [a £ {1,2,..., iV}) 






P 


bet or all parents F = |Fi, F2, . . . , Fjv j 




en 


Va 


Set of visible final state objects associated with the a}"^ parent 


x) 


ts and s« 


la 


Set of invisible final state objects associated with the a}^ parent 
Set of all visible final state objects (V = {Vi, V2, • • . ,VjVv}) 


E] 
H] 


bjec 




Set of all invisible final state objects (I = {Ii,l2, . . . ,1nx}) 


bO 


o 


iV = \v\ 

Nv = \V\ 

Ni = \I\ 

indices 


Number of parents assumed for the interpretation being applied 
Total number of visible final state objects 
Total number of invisible final state objects 

f For notational purposes, indices are used interchangably with the 

1 n n XTi pt; r»"F i"Vi a T^fi"ri"ipl pt; i" Vi p'\7" i H PTii"i PVit pypi ttit^I p* " pti H "^y^mi " p tp I 

1 llCLlilCO KjL LllC UCLi- Ijlv-lCo uliCV HJ-CllLlly- 1- VJi. CACllllL/lC r a CLllVj. ' a r, di. C 1 

1 equivalent; "i G V" and "i G {1, 2, . . . , A'"v}" are equivalent; "a £ V" \ 
I, and "a £ {1, 2, . . . , A''}" are equivalent, etc. J 


E 


a 


-^i = (J^i,PiT,Piz) 


1+3 momentum components of the i^^ final state visible object (i £ V) 




moment 




Hypothesized 1+3 momentum components of the i"^ final state invisible {i G I) 






Sum of 1+3 momentum components of visible objects belonging to parent Pa 




CO 

+ 


— 

— Z^i^x^ 


Sum of 1+3 momentum components of invisible objects belonging to parent Pa 


(571 






Total 1+3 momentum components of the 'UVM' set 


(51 1 




Pt 


Missing transverse momentum two vector (magnitude j^j, = p^) 


(51 1 




Ma = Mf^ 


Mass 01 tne a parent ya G r J 






A/f M.r 

IVli = iViYi 


IVIass of the i, visible final state object {% G V) 






IVIi = iVli. 


Hypothesized mass of the z invisible {% G ^) 




.2 


r. \ T\ ,f \ ^ r- T I 

= S Mi \ l ^ ±a} 


Set of hypothesised masses of the invisibles associated with parent Fa 




quar 


~ 1 1 ~ 
M = Ua Ma 


bet or the hypothesised masses or all invisibles 




derived 


A A 

Ma 


Hypothesized 1+3 dim. invariant mass of the composite parent particle Fa 


(62l 




1+3 dim. invariant mass of the visibles in Va 


(|o3p 




Ayr 


1+3 dim. invariant mass of the invisibles in Xa 


(54 1 




T/ 


o-speea or tne t visiDie [t v) 






V- 


nypuLiifcjbiztJu. o-bpct;ci ui Liie t iiiVibiLUc 1 1 t -i- ) 






y = /y. 1 i p T \ 

— J • i ^ ^ -^a 1 


Hypothesised 3-speeds of the invisibles associated with parent F^ 


([sol 




5 = Ua 


^pt OT ri\mnt nPQicpH qt^ppHq c\\ all trip im/'iQi rilpQ 
kJCl) VJl 11 V Mw lillColoC^J. O oL'CCU.o \jL cill Lllc; lllVlolUlCo 


(4-Q 1 




IV^ — \ M \ 
J-V^a — Z^ieia [^*-'»J 


Lllli (Ji Llic llidrObcn (Ji LliUbe illVlbiUieb cLbbUClcLLctl WlLll pcLItJllL K a 


1 »0 1 






Set of all 'invisible particle mass sum parameters' 


(941 






Largest hypothesised 3-speed of any invisible associated with parent Fq, 


/ 1 1 7 i 






Set of all 'maximum invisible 3-speed parameters' 


(1181 


1+2 d 




1+2 dim. projected energy-momentum vector for the i^^ visible 

Hypothesized 1+2 dim. projected energy- momentum vector for the i^^ invisible 


Sec. Ill 



TABLE II: Notation used in the description of events. 
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Parent I 



Qi' 



Parent P, 



Parent P3 



PS 



Parent Pi 

111 



Q5' 



FIG. 6: This figure illustrates the notation used to label 
physics objects and their assignments to parent hypothe- 
ses. The figure shows a hypothesis in which six {Ny = 6) 
visible physics objects V = {Vi, V2, V3, V4, V5, Ve} and five 
{Nx = 5) invisible physics objects I = {Ii, I2, 13, 14, 15} have 
been assigned to three (TV = 3) parents V — {Pi,P2,P3} 
according to the assignments Vi = {V2,V4}, V2 = {V3}, 
V3 = {Vi,V5, Ve}, Ii = {Ii}, I2 = {12,13} and I3 = {UM- 
The number of visible physics objects assigned to each parent 
in turn are therefore |Vi| = 2, IV2I = 1 and IV3I — 3 and the 
number of invisible physics objects assigned to each parent in 
turn are — 1, II2I = 2 and II3I — 2. 



FIG. 8: This figure is provided for the benefit of readers 
unable to imagine an even simpler version of Figure [6] than 
was shown in Figure [7] (Readers finding this figure helpful 
are advised to seek gainful employment in some other field.) 
The figure shows a hypothesis in which two (A^v = 2) visible 
physics objects V = {Vi,V2} and two (Nx = 2) invisible 
physics objects I = {Ii, I2} have been assigned to one (A'^ = 1) 
parent V — {Pi} according to the assignments Vi — {Vi, V2} 
and Ii = {Ii,l2}. For completeness we note |Vi| = = 2. 
An explicit physics example corresponding to this figure is 
discussed in Section TXI Al 



Parent Pi 



Parent P2 



> > 



FIG. 7: This figure is provided for the benefit of readers un- 
able to imagine a simpler version of Figure[6] (Readers finding 
this figure helpful need not admit this to close friends, rela- 
tives or colleagues.) The figure shows a hypothesis in which 
four (A^v = 4) visible physics objects V = {Vi, V2, V3, V4} 
and two (A^i = 2) invisible physics objects I = {Ii,l2} have 
been assigned to two (A'^ = 2) parents V = {Pi,P2} according 
to the assignments Vi = {Vi,V3}, V2 = {V2,V4}, Xi = {Ii} 
and T2 = {12}. The number of visible physics objects assigned 
to each parent in turn are therefore |Vi| —2, IV2I =2 and the 
number of invisible physics objects assigned to each parent 
in turn are \Xi\ — 1, II2I ~ 1. An explicit physics example 
corresponding to this figure is discussed in Section pCIB[ 



it is useful to have notation that can refer to composite 
quantities, e.g. the total four momentum posessed by the 
visible daughters of Pa , or the total invariant mass of that 
collection of visible daughters. Accordingly, as illustrated 
in Figures |6) [7] and |8j we denote by the total (1-1-3)- 
momentum of the visible daughters of parent Va- 



or in components 



E 



a, PaT, Paj 



PaT 



Paz 



E 



(53) 



(54) 



(55) 



E V^^- 



-pIt+pI 



- Ef -E 



h+pl 



V, 



(56a) 



(56b) 



where the former (latter) expression for will be rele- 
vant later on for T (V) transverse projections since it is 
written in a form which depends explicitly on the masses 
(speeds) of the visible particles. 

Similarly, we denote the total hypothesized (H-3)- 
momentum of the invisible daughters of parent Pa by 



Q^= (Ea,qaT,qa.) = E^^ 
or in components 



i£la 



E 9^^' 
E 



E„ 



2 

iz ' 



ieia 



V,, 



(57) 

(58) 
(59) 

(60a) 
(60b) 



where again the former (latter) expression for E^ will 
be relevant for T (V) transverse projections since it is 
written in a form which depends explicitly on the masses 
(speeds) of the invisible particles. 

As already indicated in eqs. ( 53|[60 l, we shall use bold- 
face script to label "composite" momenta. Each parent 
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is thus also treated as a composite particle, which has 
(1+3) momentum 

K + Q^a (61) 

with (a priori unknown) (1+3) dim. invariant mass 



Ma = VffM- (Pa + QaT (Pa + Qa) 



(62) 



The important distinction between the bold-face nota- 
tion for composite momenta and the ordinary notation 
for the momenta of individual particles is pictorially il- 
lustrated in Figures [6l 7] and [8] 



Note that in (157) (whose tilde is necessary to dis- 
tinguish it from the energy Eq of the visible composite 
daughter of parent a) might legitimately be termed the 
missing energy^^ of the parent Pq. We also introduce 
masses for the respective composite daughter objects as 
follows 



Ma(Qr) 



E2 



PaT Paz! 



n2 



(63) 
(64) 



where again a tilde refers to the invisible object. Note 
that the invariant mass (64) of any composite invisible 
daughter is not "constant" or a fixed function of mea- 
sured momenta. It depends on the hypothesized invisible 
momenta Q'' and so is part of the event hypothesis. 



C. Notation for "early" and "late" transverse 
projections 

When forming transverse kinematic variables corre- 
sponding to composite parent or daughter objects, one 
needs to construct the transverse 1+2 dim. analogues of 
(53) and (57). In doing so, one inevitably has to face the 
issue discussed in Section |V] — whether the agglomera- 
tion of individual particles into a composite object is done 
before or after projecting into the transverse plane. As 
we already saw in Section|V) the two outcomes are gener- 
ally quite different, since the composite object is usually 
massive. This is why we shall need to develop some ad- 
ditional notation to help us keep track of the order in 
which those operations are performed. Correspondingly, 
for the remainder of this paper we shall adopt the follow- 
ing principle: in forming transverse quantities for com- 
posite objects, the order in which the various operations 
of agglomeration and projection are taken will be speci- 
fied by the order (from left to right) of the corresponding 
subscript indices. 

Let us illustrate this principle with a few relevant ex- 
amples. The "late-projected" (or "early-partitioned") 



version of the composite visible momentum ( 53 ) is de- 
noted by p"rp 



PaT = (ear, Par) 



(65) 



while the alternative "early-projected" (or "late- 
partitioned" ) version is denoted by p^^ : 



PTa = (eTa,PTa). 



(66) 



We remind the reader that the generic index "T" in ( 65 ) 
and (66) stands for either "T", "V" or "o", as discussed 



in Section rrni 

We already saw in Section |v] (eqs. ( 37p9 )) that the 
space-like components of (|65| and (|66]) are equivalent for 
any choice of "T" : 



PaT = PTa=^ PiT , 



(67) 



but the time-like components BaT and Bxa are gener- 
ally different. For example, in the case of T = T, the 
late-projected (early-partitioned) transverse energy e^T 
is given by 



y^l + Plr-^K-pl. 



\ 



^+pfT+Pl 




(68a) 



,(68b) 



while the early-projected (late-partitioned) transverse 
energy era is given by 



(69) 



era = ^'^ = E V^^'+PtT 



In the case of T = V projections, the corresponding trans- 
verse energies are given by 



e„v = 



PaT 



VPaT + PI: 



E„ 



E\ - ViT 
= 2^ IT ■ 

Finally, for T = the two transverse energies are 



Ep»- 



(70) 
(71) 

(72) 
(73) 



The same conventions apply to the transverse projec- 
tions of the composite momentum of a collection of invis- 
ible daughter particles: the "late-projected" (or "early- 
partitioned") version of the composite invisible momen- 
tum (57) is denoted by q"j. 



Really the missing energy rather than the missing momentum}. 



(74) 
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while the alternative "early-projected" (or "late- 
partitioned") version is denoted by q^^,: 



QTa = (,eTa,qTaJ 



(75) 



Again, the space-like components of ( 74 ) and ( 75 ) are the 
same: 



(76) 



but the time-like components are not. Altogether, there 
are 6 different possibilities: 




era 



?2 



E 



E QiT 



E 



(77a) 



(77b) 



(78) 



(79a) 



(79b) 



(80) 



(81) 



(82) 



9 N = 2. This assumption is motivated if the new 
particles are charged under a conserved Z2 par- 
ity, like i?-parity in supersymmetry or KK-parity 
in UED. However, other discrete symmetries are 
also possible, e.g. Z3 [101 [H] and higher 
which could in principle allow for N > 2. Even 
in models with a Z2 parity one could still consider 
the production of any even number of parents, e.g. 
N ^4, N = 6, etc. 

• \Ta\ = 1 for all a. In the conventional models with 
conserved Z2 parity, this assumption implies that 
the decay of each parent generates one and only 
one massive invisible particle, excluding the pos- 
sibility of any neutrinos appearing among the in- 
visible particles. However, this assumption is not 
guaranteed — even in the conventional SUSY mod- 
els with conserved i?-parity, SM neutrinos can eas- 
ily appear among the decay products of charginos, 
sleptons, T4^-bosons, heavy flavor quarks (especially 
top), taus, etc. Furthermore, i?-parity conservation 
only guarantees that a given SUSY parent must 
decay into an odd (not necessarily 1) number of 
SUSY particles. Finally, a Z3 symmetry could al- 
low two massive invisible particles per parent, see 

e.g. 01 sg. 

Because of all these caveats, we prefer to keep our discus- 
sion as general as possible, and first define our invariant 
mass variables below in Sec. |VII| for any TV and Nx, be- 
fore specializing to = 1 and A = 2 for illustration 
purposes only. 

One might ask whether the methods proposed here can 
be usefully applied to events with "the wrong" value of 
A^. The answer to this question is "yes", and we shall 



demonstrate this explicitly below in Section XI B (see 
in particular Figure 12 ) where we shall apply A^ = 1 



variables in an example where not one, but two parents 
were produced in the hard scatter. That study will show 
that one can sometimes obtain useful information from 
variables with "the wrong" value of the parent number 
N. 



In general, our principle of index ordering will extend 
to any transverse invariant mass or transverse energy 
variable. For example, in analogy to (68 73) and (77 



82 1 there will be six different versions of the transverse 
masses of the composite parent particles and they will be 
denoted by MaT or Mra, with T e {T, V, o}. 



D. Comments on the characterization framework 

Note that we do not place any a priori restrictions on 
the values of A^, Nx or on the way invisible particles 
are partitioned into the subsets la- In contrast, many 
studies on supersymmetry (SUSY) or Universal Extra 
Dimensions (UED) in the hadron collider literature are 
predicated on the following assumptions: 



E. Choosing the partitioning 

In conclusion of this section, one more comment re- 
garding the partitioning is in order. One may wonder how 
one should decide whether a given visible particle should 
be counted among the set of visible daughters or whether 
it should be included in the "Upstream visible momen- 
tum" category. The answer to this question depends on 
the particular case at hand. There are simple cases of 
final states where the outgoing particles can be unam- 
biguously associated with the particle sets Va that match 
the expected decay products of an assumed parent. For 
example, a high px, isolated reconstructed lepton is un- 
likely to have come from the typical sources of UVM such 
as initial state radiation (ISR), multiple parton-parton 
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interactions (MPI), multiple hadron-hadron interactions 
(pileup) etc., and can probably be safely counted as a 
visible daughter. On the other hand, there are also cases 
(typically involving jets of hadrons) where the correct 
partitioning is not obvious at all. In such cases, one pos- 
sible approach is to consider all possible partitions, see 
e.g. [3. Another possible approach would be to devise a 
certain set of cuts, using the generic differences between 
the kinematics of ISR jets and jets from heavy parent 
decays [28l I46ti48] . Examples of choices for particular 
physical examples can be found in Section |X) 



VII. THE MASS-BOUND VARIABLES 
A. Guiding principles 

The guiding principle we employ for creating useful 
hadron-collider event variables, is that: we should place 
the best possible bounds on any Lorentz invariants of in- 
terest, such as parent masses or the center- of-mass en- 
ergy s^/^, in any cases where it is not possible to deter- 
mine the actual values of those Lorentz invariants due to 
incomplete event information. Such incomplete informa- 
tion could take the form of lack of knowledge of the lon- 
gitudinal momentum of the primary collision, or lack of 
knowledge of the 4-momenta of individual invisible par- 
ticles, or lack of knowledge of the number of invisible 
particles which were present, etc. 

We contrast this principle with the alternative ap- 
proach that is used to motivate event variables without 
any explicit regard to whether they have an interpreta- 
tion as an optimal bound of a Lorentz invariant. This 
alternative approach tends to recommend the use of vari- 
ables that are somewhat ad-hoc, but by construction pos- 
sess useful invariances (such as invariance under longitu- 
dinal boosts) which are designed to remove sensitivity to 
quantities that are unknown. One example of this lat- 
ter class of variables, which are usually considered to be 
simply "made up" without reference to our guiding prin- 
ciple, would include the missing transverse momentum 
firp (aready seen in (51 1 ) obtained by adding all trans- 
verse visible momenta vectorially. Another would be the 
so called Ht variable'^^ which is defined as the scalar sum 



of the transverse momenta of some class of visible objects 
(typically jets) in the event: 

hT = J2P^T- (83) 
Another example is the sum of these two variables: 



(84) 



a quantity which can be traced back to the original liter- 
ature [49J and has become known as an "effective mass" , 
even though it is not a mass.^'^ 

The main disadvantage of variables like Ht and TOcff , 
is that they do not utilize all the information available; 
for example they are completely insensitive to all angles 
in the transverse plane. This is why here we would like 
to construct a more optimal class of variables, to wit, 
those which bound the invariants of interest. These too 
must be invariant under global longitudinal boosts since 
a bound cannot depend on unknown quantities. However 
by explict construction we can ensure that they also make 
best use of any available kinematic information. 



B. Construction of mass-bounding variables 

We are now ready to define the general procedure that 
can be used to construct the mass-bound variables. In 
fact, we shall describe a broad class of such variables, 
where each individual variable M^indices} will be labelled 
by a certain set of indices {indices} indicative of the way 
the particular variable was constructed, namely: 

• Since we are targeting the general event topology 
of Fig. [5j where we imagine the inclusive produc- 
tion of N parents, each one of our variables will 
necessarily carry a corresponding index N. In the 
process of constructing such a variable, we will 
have to partition (and then agglomerate) the ob- 
served visible particles in the event into N groups 
Va, {a = 1,2, ... ,N), as already explained in Sec- 
tion [VI^ We will then form the 1-1-3 dimensional 
invariant mass of each parent 



^■^ Note that the definition of Ht in the literature is not well stan- 
dardized. Indeed even one LHC experiment has managed to 
define it in three different and inequivalent ways in the space of 
just a few years, and sometimes even inconsistently in a single 
document (see Section 2 of U for further details). The defini- 
tion we adopt in equation \83\ is the definition which appears, 
at present, to be the most widely used in the literature. We 
note that a conceivable consequence of this paper might be that 
purists will in the future settle on a definition in which is 
defined as a sum of transverse energies st instead of transverse 
momenta, whereby three different variants would be possible: 
hy, hv and ho (though these three definitions will be almost 
equivalent under most practical experimental conditions, where 



Jg^,(P^ + Qj^)(P- + Q;i), 



(85) 



which is constructed out of the 1+3 momenta 
and of the respective composite daughter par- 
ticles (see Section VI B ) . 



the visible particles are approximately massless). 

In keeping with our conventions from Section|ll] we use lowercase 

letters for both hx and rricS, since they are not H-3 dimensional 

quantities. 
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Type of 
variables 


Operations 


Notation 


First 


Second 


Third 


Unprojected 


Partitioning 


Minimization 




Mjv 


Early partitioned 
(late projected) 

MjVT 


Partitioning 
Partitioning 
Partitioning 


T = T projection 
T — y projection 
T — o projection 


Minimization 
Minimization 
Minimization 


Matt 
Mnv 

Mno 


Late partitioned 
(early projected) 

Mtn 


T = T projection 
T = y projection 
T = o projection 


Partitioning 
Partitioning 
Partitioning 


Minimization 
Minimization 
Minimization 


Mtn 

MvN 
MoN 



TABLE IIL Method of constructing the mass-bound variables and corresponding notation. The Table lists the sequence of 
operations performed in the calculation of each variable. "Partitioning" refers to the operations discussed in Sec. |VI A| and 
I VI B| of partitioning the final state particles into daughter sets and then adding the momenta in each set to form corresponding 
composite daughter particles. "Minimization" implies minimizing the largest (suitably projected) parent mass with respect to 
(the relevant components of) the missing momenta of all invisible particles; while the remaining operations involve the different 
types of transverse projections defined and discussed in Section |III[ 



• Optionally, instead of the 1+3 dimensional parent 



mass ( 85 ) , we may choose to consider the corre- 



sponding early-partitioned (late-projected) trans- 
verse mass 



Ma7 



5o^(P^T + qaT)(pfT + qaT): (86) 



or the late-partitioned (early-projected) transverse 
mass 



M 



Ta 



M(PTa + q?a)(PTa + qTa). (87) 



/9 



where p"^, p^^, q"^ and qf.^ are the 1+2 dimen- 
sional momentum vectors defined in (|65|), ([66|), ([74|) 



and ( 75 ) , correspondingly, and the index T t akes 

14 



values in {T, V, o}, as explained in Section 



III 



We should point out that the projection specification T £ 
{T, V, o} refers to operations on the visible particles. One should 
keep in mind that the visible and the invisible composite parti- 
cles are a priori independent and so could, in principle, be treated 
differently, both in terms of the order of the operations, as well 
as regarding the type of transverse projections. For example, 
consider the M]vt class of variables, where one first forms com- 
posite visible particles and transversifies later. In principle, for 
the invisible particles, one could perform those operations in the 
opposite order and instead of ||86|l consider 



instead. Furthermore, one could choose a different type of trans- 
verse projection for the invisibles than for the visible sector, e.g. 



(Pal 



/)(pfT 



and so on. One might therefore wonder whether projected vari- 
ables need to carry additional indices indicating how the invisible 
sector is being handled. In the following, for simplicity we shall 
assume that the invisible particles are always projected in exactly 



• The last step is to consider the largest hypothesized 
parent mass (max[A^a], max[A^QT] or maxfA^^a] 
as appropriate) and minimize it over all possible 
values of the unknown invisible momenta consistent 
with the constraints. This minimization is always a 
well-defined, unambiguous operation, which yields 
a unique numerical answer |50) . which we shall de- 
note as 



M 



N 



M. 



NT 



TN 



mm 



mm 



mm 



max \M.c\ 



max [A^qt] 



-\M 



Ta\ 



(89) 



(90) 



as indicated in Table lllll The minimization over the 
unknown parameter is performed in order to guar- 
entee that the resultant variable cannot be larger 
than the mass of the heaviest parent, resulting in 
an event-by-event lower bound on the mass of the 
heaviest parent. 

These are the basic steps, leading to the variables dis- 



played in Table III This basic set of variables will be 
further extended in Section |VIII| below, by considering 
a second level of projections within the transverse plane. 
For the remainder of this section, however, we shall stick 
to the basic procedures above and focus on the simplest 



the same way as the corresponding visible particles, so that the 
transversification indices uniquely describe the transverse pro- 
jections of both visible and invisible daughters. Those readers 
who are curious about the remaining cases (when the visibles 
and the invisibles are projected differently) can easily infer the 
corresponding results from the formulas given below. 
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classes of variables displayed in Table namely the 
"unprojected" Mm and the "singly projected" Mnt and 
Mtn variables. 



C. The variables: Mjv, Mnt and Mtn 

In this subsection we provide analytic formulas (where 
available) for calculating each of the basic mass-bound 
variables from Table [IIT| on an event-by-event basis. 



1. The usual ( "unprojected" ) invariant mass: Mm 

Here we work directly with the usual (1-1-3)- 
dimensional invariant masses A^q of the parent particles 



Ml{Va,Qa,K) = {'Pa+Claf 



(91a) 



= (Ea + E„)' - (paT + qar)" - (Pa. + Qa^)" • (91b) 

The unprojected invariant mass variable Mjv is defined 
by the right hand side of 



Mjv(]yi) = niin 



max[A^a(Pa,Qa,Ma)] 



(92) 



where the minimization needs to be performed over 37Vx 
degrees of freedom {qiT and qiz for i = 1, 2, ... , Ni), 



subject to the two scalar constraints (51 1 supplied by 



transverse momentum conservation. The invisible parti- 
cle momenta qiT and qiz are fixed by the minimization 
and does not depend on them. 

Note that we have emphasized in the left hand side of 
( 92 1 that Mn turns out not to be a function of the Ni 



individual invisible mass hypotheses Mi in /i = IJaAa, 
but instead turns out (sec proof in Section IX A I to be a 
function of the set 



(93) 



containing the N "invisible mass-sum parameters, IV^o" 
defined by 



(94) 



These mass parameters are simple arithmetic sums of the 
hypothesized masses of the individual invisible particles 
associated with any given parent Pq. 

Notice the simplification in going from the individual 
parent masses M.a to the variable Mm- The individ- 
ual parent masses AAa collectively depend on all invisi- 
ble particle masses Mi, (a total of Ni parameters), while 
the invariant mass variable Mjq defined in ( [92| only de- 
pends on the N summed-invisible-mass parameters IV^ai 
(a = 1,2,..., TV), which are simply related to the indi- 
vidual particle masses Mi via ( 94 ) . In the most common 



cases of A'' = 1 or iV = 2, we will therefore have to 
deal with only one or two unknown invisible mass-sum 
parameters. A similar reduction in complexity will be 
found when we consider the V projected variables, but 
there the mass bound will end up depending on a speed- 
related parameter for each parent. We see that from now 
on the index N can be interpreted not only as the number 
of parents, but also as the number of relevant indepen- 
dent mass inputs characterizing the invisible sector. 

The preceding discussion is best illustrated with a spe- 
cific example. Let us consider the simplest case of TV = 1. 
The minimization of the corresponding variable Mi with 
respect to qiT and qiz is straightforward. One finds that 
the minimum is located at [301 



qiT = Pt 



Mi_ 
M 



and its value (see [H]) is given by 



(95) 
(96) 



m| (97) 



in which, to save space, we have slightly abused our no- 
tation by writing Mf(]yii) in place of Ml{{^i)) — 
a convention we will adopt throughout this document 
wherever N — 1. We remind the reader that Mi is the 
measured (l-|-3)-mass of the (single) visible composite 
daughter (see also eq. (63)) 



Ml = J^l - - p2 



(98) 



while lyii is the only invisible mass parameter needed^^ 
defined in (Igll) 



Nx 

1^1 EE ^Af, 

i=l 



(99) 



In Ref. [46], the quantity Mi{lfii) defined in ([97]) was 
labelled "\/ii,;„ ^ 



Ml (1^1 ) 



-(sub) 



(100) 



since it provides a lower bound on the parton-level center- 
of-mass energy of the parent subsystem Vi © Ii, not 



1^ Note the analogy between fj, and l^Ii. frp measures the total 
transverse momentum of the whole collection of missing particles. 
Similarly, 1^1 measures the total mass of the whole collection of 
missing particles. Both pj, and ifii are given by simple sums 
of the corresponding quan titie s QiT an d Mi of the individual 
missing particles, compare (|51|l and \99\. 
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counting the uninteresting upstream visible momentum 
U'^ . In the special case of a vanishing upstream momen- 
tum {ut — 0), Afi(]yii) reduces to the global variable 
Vsmin from [50]: 



lim Mi(]yii) 



(101) 



We will not consider the next simplest example {M2) 
until Section XH as simple analytic (as opposed to nu- 



merical or iterative) formulae for it are only known to 
exist in certain special cases [5T], such as when l^I^ — 
1^2 — Ml = M2 — 0, or when the upstream visible mo- 
mentum Ut is either zero or (anti-)parallel to the missing 
transverse momentum . 



2. The early partitioned, T -projected invariant mass: Mnt 

Here the momenta and of the composite par- 
ticles are first formed in 1-1-3 dimensions, as in (53) and 



(571, then afterwards are projected on the transverse 



plane according to the mass-preserving T method de- 
fined in cq. 



of Sec. Ill A This results in transverse 



masses of the parents given by 

A^aT(PaT,qaT,Aa) ^ (PaT+Qar)' (102a) 
= {ea-T + ea-T? - {^aT + ^arf , (102b) 

where the transverse momenta and (\aT are given by 
(54) and (58), while the transverse energies GaT and e^T 



are given by ( 68 ) and ( 77 ) . 

Then the "early partitioned, T-projected" variable 
MjvT is defined in a manner very similar to ( 92 ) 



Mwt(¥I) 



mm 



max[7WaT(PaT:qaT>Aia)] 



(103) 

Just like Mjv, this variable also depends only^^ on the 
N summed-invisible-mass parameters lyl^j within 1^1 as 
opposed to the iVx individual invisible masses Mi within 
/i. Eq. ( 103 ) again represents a constrained minimization 
problem for the 3A'x variables qiT and qiz- Note that in 



^® At this point readers who are familiar with the Cambridge 

variable [71[8| have probably recognized that for the special case 
of A'^ = 2, the MjvT variable |103| recovers the Cambridge mT2- 
Note that the original literature ^Sj on the Cambridge mx2 vari- 
able also defined more general variables rriTX, e.g. mx3, rriTA, 
etc. However, we caution readers to make the distinction between 
the index "A" in M]vt, which refers to the number of hypothe- 
sized parents, and the index "A" in the Cambridge rrixx i which 
stood for the total number of invisible particles (in this paper 
denoted by Ax). For example, the index "2" in the Cambridge 
mT2 notation implies the presence of exactly two invisible parti- 
cles, the number of parents already being implicitly assumed to 
be two. In contrast, the variable M2T defined in ( |103[ l does not 
imply any particular number of invisible particles, and in this 
sense is equivalent to the whole class of rnj^x ior any X. 



spite of its "transverse" index, Mnj- still depends on the 
longitudinal momenta Qiz through the transverse energy 
BaT, see (|77b|). 



In order to gain some intuition, let us again consider 
the simplest case of = 1. The minimization of ( 103[ ) is 
once again straightforward and the minimum is found at 



qtz = qiz 



(104) 
(105) 



with an arbitrary choice of qi2. This leads to 



Ut ■ 



Comparing (106) to (97), we see that 



MiT = Ml. 



(106) 



(107) 



This is in fact a special case of the more general mathe- 
matical identity 



M 



NT 



N, 



(108) 



for which a proof is provided in the appendix — see equa- 
tion ( B49 ) . This identity reveals that "transverse" quan- 



tities do not necessarily "forget" about relative longitudi- 
nal momenta. In particular, ( 108 ) teaches us that when- 



ever the composite particles are formed before the trans- 
verse projection, the information about the relative lon- 
gitudinal momenta is retained, and the result is the same 
as if everything was done in 1-1-3 dimensions throughout. 
As a result, Mnj- automatically inherits all the advan- 
tages and disadvantages of its 1-1-3 cousin Mjy. 



3. The late partitioned, T -projected invariant mass: Mtn 

This is the first example of an "early projected" , "late 
partitioned" variable. We follow the procedure of the pre- 
vious subsection I VII C 2] only this time we switch the or- 
der of the operations, and we first T-project the momen- 
tum of each individual particle on the transverse plane, 
before forming composite particles. The transverse in- 
variant mass of each composite parent is then given by 



XTa(PTa>qTa:Ma) = (PTa + qTa) 

= (era + BTq)^ - (PaT + ^aT. 



(109a) 
2, (109b) 



with and qaT still given by ( 54 ) and ( 58 ) , while the 



composite transverse energies 6-^^ and e-^a are given by 
(69) and (l78|), correspondingly. Notice that these ex- 



pressions do not contain the longitudinal momenta piz 
and Qiz- This is in contrast to the "early partitioned" 
case represented by ( 68 ) and ( 77 ) , where the longitudi- 
nal momenta appear explicitly. The comparison between 
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(681 and (77) on the one hand, and (69) and (78) on the 
other, nicely illustrates the main point of Section |V] — 
that by adding the momenta before the projection, one 
retains sensitivity to the relative longitudinal momenta. 
Conversely, when the operations are performed in reverse 
order and the transverse projection is done first, the lon- 
gitudinal momenta completely drop out of the game. 

Now we are ready to apply the usual definition and 
obtain 

MTAr(]yi) = min max [AlTa(PTa: QTa- Aa)] 

(110) 

Let us again investigate the simplest case of iV = 1. 
With the help of the transverse momentum conservation 



constraint (|51|), eq. (|110| reduces to 



mm 

E <iiT=fT 



\i=\ 




:i=l 



The minimum is once again found at ( 104 ) and we get 




(111) 

As expected, this result differs from (1061, although the 
two formulas follow a similar pattern. The difference 
is only in the term corresponding to the visible sector, 
where the transverse energy of the composite visible par- 
ticle is computed differently, compare ( 68a ) and ( 69 1 . 

An interesting result emerges if we consider the fur- 
ther simplification that all visible particles are massless, 
i.e. Mi = 0,V«. This, in fact, is a very good approxima- 
tion for the leptons and quarks/gluons of the SM, whose 
masses can be safely neglected. Setting Mi = in (111) 
and using ( 83 ) , we get 



lim M^i(]yii) = hr 



mi 



(112) 



This result is quite interesting. It allows us to reinterpret 
the usual hr variable in terms of a bona fide invariant 
mass variable like Mti, properly accounting for the ef- 
fects of upstream visible momentum ut and the total 
mass of the invisible particles present in the event. 
We shall return to this point in the next Section 



IX 



Another interesting result follows from eq. ( 112 ) in the 
special case when we set IVli = 0. Using (84), we get 



Jni^A//2,(D/li = 0) = (^IiT+^^y -u\ = mis - 4, 

providing a connection between the "effective mass" TOoff 
and Mti(0). 



4- The late partitioned, V -projected mass: M\/n 

This is the second example of an "early projected" vari- 
able, only this time we use the speed-preserving V pro- 
jection described in Section |IIIB Correspondingly, the 
individual visible (invisible) particles will be character- 
ized by their 3-speeds Vi {Vi) instead of their masses Mi 
(Mi) and so we remind the reader of the notation intro- 
duced in (|49l) and (1501). 



The 1-1-2 momentum vectors of the individual particles 



after the V projection are obtained from (17) 



Pfv = (eiv,Kv) = y-^,PtTj , (114) 
Qiw = (ejv,gjv) = (^^,giT^ • (115) 
Then we form composite particles with V projected 1-1-2 



momenta Py^ and qy^ given by (66) and (75), respec- 
tively. 

The transverse parent masses are now formed in terms 
of Py^ and qyjj as follows 

^laiPva^<iva^Va) = {pva+^lvaf (116a) 
EE (eya + eya)^ - {paT + ^axf, (116b) 

where the transverse energies eya and eya are specified 
by ( 71 ) and ( 80 ) and the transverse momenta PaT and 



QaT are given by ( 54 ) and ( 58 ) 



This is a convenient place to introduce another two 
small pieces of notation. Firstly we will need to define 
a "maximum invisible velocity parameter" Va for each 
parent Pa according to 



Va = max 

ieia 



(117) 



Then we would like to denote by V" the set of all the 
above velocity parameters, i.e. 

y={ya|aeP}. (118) 

Now we are in a position to state (see proof in Sec- 



tion 



IX A) that the only dependence of the "late parti- 



tioned", V-projected mass variable My^r on the velocity 
parameters of the invisible particles is through V, i.e.: 



A/viv(y) 



mm 



max[A^ya(Pva,qva,Wa)] • (119) 



Once again, it is instructive to consider the special case 



Contrast with the definitioi i of iVl^ in equation l |94| and the 
definition of ly^ in equation |93|. 
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of = 1. With the help of eq. (|TT9| becomes 

2 




The minimization selects the invisible particle with the 
largest speed, whose transverse momentum becomes ^j,, 
while all other invisible particles have qiT = 0. This 
configuration leads to the final answer 

^ 2 



(120) 



When we make the approximation that all visible par- 
ticles are massless {Vi = 1), we again obtain a relation to 



fT 



(121) 



which is the analogue of ( 112 ) for the case of V transverse 



projections. But note that unlike (112), here the un- 



known parameter characterizing the invisible sector is the 
maximum speed parameter V'j^ instead of the summed- 
invisible-mass parameter iVl^. 

Finally, if in addition we also assume that all invisible 
particles are massless as well, then 



V, = 1, Vz 



^1 = 1, 



so that 



'cff 



^ ^ ^ ^(122) 

which is the analogue of (1131. The fact that ( |113| ) and 
(122 1 arc the same should not come as a surprise: recall 



from Sec. IV B that the two transverse projections T and 
V are equivalent in the massless limit. 



5. The early partitioned, y -projected mass: Mnv 

Here we follow a procedure analogous to that of 
Sec. VII C 2 where the composite momenta and 



are first formed in 1+3 dimensions, before being pro- 
jected on the transverse plane, only this time we use the 
V projection for this purpose: 



= (eaV,Pav), (123) 



qav = (eav,qav). (124) 



The transverse parent masses are now formed in terms 
of p"v and q"v as usual 



(125a) 



(eav + ea 



(PaT + qaT)'(125b) 



Here the composite transverse momenta p^v and q^v are 
still given by ( 54 ) and ( 58 1 , while the transverse energies 



Gav and Bav are given by ( 70 ) and ( 79 ) , correspondingly. 



Then the early-partitioned, V-projected variable is de- 
fined as usual: 

Mnv i'f)= min max [M av (Pav > Qav > )] 

J2<liT=pT L " 

(126) 

Once again, let us specify this to the case of iV — 1. 
Using (51), we get 



Ml, 



mm 



[eiv + eiv) 



(127a) 



= eiv 



min [eiv] ) - (127b) 



The minimization is performed over the 3Nx variables 
QiT and Qiz, i ^ 1,2, . . . , Nx and the result is 



Miv(yi) = e 



Pit 



El 



(128a) 



- u^,(128b) 



which is similar, but not equivalent to ( 120 ) 



6. The late partitioned, o -projected mass: Mo 



Here we follow the procedure of Sees. |VII C 3| and 
VII C 4[ only this time we use the o transverse projec- 
tion from Sec. |III Cl One first forms the 1-1-2 momenta of 
the individual particles 

Pfo = [eiojio) = [piT, Pit) , (129) 

qfo = {eio,qto) = {qtT,qiT) , (130) 

then the composite momenta 

Poa = (eoa,Poa) = ^ E ' ^^^^^ 

\ieVa ieVa / 

Qoa = (eoQ,qoa) = E E ) • 

The transverse parent masses are now formed in terms 
of Po£, and q"^ as usual: 



-^oa(Poa.qoa) = (Poa + qoa)' 



(133a) 



(eoa+eoa)^-(PaT + qaT)^ (133b) 
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and the "late partitioned", o-projected mass variable 
MoN is defined as before: 



MoN= min max [7Woa(Poa> Qoa)] 



(134) 



Notice that the Mon variables do not depend on any un- 
known parameters related to the invisible sector (i.e. we 
need no "0" where previously we needed an l^A or a ^) 
and so can be uniquely computed in terms of the mea- 
sured momenta of the visible particles and the missing 
transverse momentum alone. 



Specializing (134 1 to the simplest case of = 1, we 



get 



mm 





mm 

12 ?iT=pT 



mm 



'Nx 

.1=1 

'Nx 

E 

i=l 



qiT 




The minimization over the 2Nx variables qiT is straight- 
forward and we obtain several equivalent expressions for 
the answer 

2 

M'i, = I > ^1D,T + J^„ 1 -4 (135a) 

(135b) 

— miff — u^. (135c) 

showing the close connection between Moi and the usual 
and rrics variables. We see that in the absence of any 
upstream visible momentum (ut = 0), the variable Moi 
itself is nothing but the effective mass rricff- However, 
these two variables differ if (as is typically the case) the 
event also has some nonzero upstream momentum ut- 
The importance of the result ( 135c I is that it teaches us 



how to properly account for the presence of UVM in such 
cases: ut should be subtracted in quadratures from rUeS 
in order to obtain the proper invariant mass variable (in 
this case Moi). Furthermore, it also reveals the phys- 
ical meaning of the widely used rUes variable (see also 
Sec. XB below): it is the minimum allowed transverse 



mass constructed out of "o"-projected momenta, for a 
semi-invisibly decaying parent, whenever that parent is 
produced exclusively with ut = (i.e. with no additional 
upstream momentum in the event). 



7. The early partitioned, o -projected mass: Mmo 

Finally, we discuss the early partitioned, o-projected 
version Mno, where the composite momenta are first 



formed in l-|-3 dimensions, then transversified via the 
"o" projection: 

J2 n ^ - (eacPao) , (136) 



Eq^ 



^ q:o = (eao,q„o), (137) 



where in light of (72) and (81 1 



(^ao ; Pao ) 



E*^ 



, E 



' E "i^T 



(138) 



(139) 



These (1+2) composite momenta are now used to form 
the corresponding transverse parent masses 

-^ao(Pao,qao) = (Pao+qao)' (140a) 
= {eao + eaof - (PaT + ^arf ■ (140b) 

Now the "early partitioned" , o-projected mass variable 
Mno is defined as before: 



Mno — niin 



max[X,o(p^o,q:o)] • (141) 



Just like its cousin Mqj^ defined in (134), M^a does not 
depend on any unknown parameters like lyl^ or 

Specifying ( 141 ) to the simplest case of = 1, we get 



— min 



(eio-l-eio) -U' 



eio + mm [eio 



eio + mm 



Nx 



4=1 



iT 



The minimization over the 2Ni variables qiT gives 

V 2 



Mi^o 



Nv 

i=i 



T + "UTl 



(142a) 



(142b) 



= 2 (/j, • {0r^ + ut) +P^\0t + ut\), (142c) 



T \yT 

providing a connection between our Afio variable and the 
usual missing transverse momentum p^. In order to see 
the physical meaning of , let us take the "no upstream 
momentum" limit ut in (142b) or (142c), resulting 
in 



lim Mf, =4^2 



(143) 
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AT, 



\ i 


Til 
^11 

1 >■ 








> 




Ut 






Transverse vector decomposition 


onto the direction 



FIG. 9: 

T|| specified by the UVM transverse momentum vector ut and 
the direction T± orthogonal to it [251 134) . All vectors shown 
are in the plane perpendicular to the beam axis. 



One can thus interpret the variable (and not just 
the as the minimum allowed "o" -projected trans- 

verse mass of a semi- invisibly decaying parent, whenever 
the parent is produced exclusively with ux = 0, i.e. with 
no additional upstream momentum in the event. How- 
ever, in situations when the parent is produced inclu- 
sively, with Ut ^ 0, the relevant variable to consider 
would be Mio as given by (142b) or (142c), which prop- 



erly accounts for the ut effect (see also Sec. XA below) 



VIII. 



ADDITIONALLY PROJECTED 
VARIABLES 



A. Momentum decompositions with respect to ut 

An additional level of projection within the plane 
transverse to the beam has been shown to be useful in 
certain circumstances [25l [34]. To orient such projec- 
tions we note that the total transverse momentum ut 
of the UVM category breaks the rotational symmetry of 
the transverse plane and selects two preferred directions 
T|| (along Ut) and Tj_ (transverse to ut), as shown in 
Fig. [9j Having projected the l-f3 momentum vectors 
onto the transverse plane as in Fig. [2j one may then ad- 
ditionally project the resulting l-|-2 transverse momen- 
tum vectors onto these special directions, as illustrated 
in Fig. [9j The corresponding momentum components 
resulting from such "double transverse" projections will 
carry a "double transverse" index: "T_L" for components 
along Tj_ and "T||" for components along Ty (see Fig. [9|. 
For example, the prp vector can be decomposed into a i± 
component f>T± 



{pt X Ut) , 



-^ut X \ pT X Ut 



(144) 



By definition, the upstream transverse momentum vector 
Ut has only a Tii component, i.e. 



UT± 



Ut, 
0. 



(146) 
(147) 



T± direction reads 

N-r 



In view of ( 51 ) and ( 147 1 , momentum conservation in the 



^ qiT± = Pt± = ^ X! Pi"^^ ■ 

i=l i=l 



(148) 



It is precisely the absence of a ut_l term in this equation 
which allows one to derive exact analytical formulas for 
the Tl doubly projected variables defined next in Section 

ivniBi 



B. Doubly-projected mass bound variables 

1. Homogeneously-doubly-projected mass bound variables 

For our purposes, the additional projections in Fig. [9] 
allow us to extend the original set of mass-bound vari- 
ables from Table |III| by consi deri ng the "doubly pro- 



jected" variables shown in TablellV^®. The benefit of such 
additionally projected varibles has been noted and dis- 
cussed in [25l[34]- For example, the shapes and the kine- 
matic endpoints of the distributions of TL-projected vari- 
ables can be independent of the value of ut- Therefore, 
such distributions can be constructed from the whole 
event sample, without any loss in statistics due to a spe- 
cific Ut selection. Furthermore, the relation (147) leads 



to significant simplifications in the analytical treatment 
of T± doubly projected variables. For example, for singly 
projected variables, the case of = 2 is untractable 
by analytical means, and (apart from some special cases 
[ST]) has to be treated numerically [351 [2S]- In contrast, 
one can derive exact analytical formulas for calculating 
N = 2, T± doubly projected mass bound variables on 
an event-per-event basis, without any need for numerical 
minimizations [3S1 [H] . 

In general, the variables in Table [TV] are independent, 
with one exception: 



M 



NoJ_ 



(149) 



Later on in Section XI B (see in particular Fig. [I3{^b)), 
we shall consider a specific example illustrating some of 
the homogeneously-doubly-projected variables from Ta- 
ble HVl 



and a T|| component ^jin 



-fT-f. 



T_L 



Ut ]ut ■ 



(145) 



To save space, Table [Tv] lists only T±_ projected variables. An 
analogous set of Ty projected variables is obtained by replacing 
the Tx projection in Table 11 V| with a Tii projection. 
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Type of 
variables 


Operations 


Notation 


First 


Second 


Third 


Fourth 


Early partitioned 
doubly projected 

MjvT± 


Partitioning 
Partitioning 
Partitioning 


T — T projection 
T = V projection 
T = o projection 


_L= T projection on T± 
_L= V projection on T± 
_L= o projection on T± 


Minimization 
Minimization 
Minimization 


A'/jVT± 

Mnv± 

MNa± 


Late partitioned, 
doubly projected 

Mt±n 


T = T projection 
T = V projection 
T — o projection 


_L= T projection on T± 
_L= V projection on T± 
A_— o projection on T± 


Partitioning 
Partitioning 
Partitioning 


Minimization 
Minimization 
Minimization 


Mt±n 

Mv±N 

Moxjv 


In-between partitioned, 
doubly projected 

Mtn± 


T = T projection 
T = V projection 
T — o projection 


Partitioning 
Partitioning 
Partitioning 


_L= T projection on Tl 
_L= V projection on T± 
_L= o projection on T± 


Minimization 
Minimization 
Minimization 


Mtjv± 

MvN± 

Mojvx 



TABLE IV: An extended version of Table |III| containing the additional variables found by including the option of a T± 
projection shown in Fig. [91 An analogous set of variables is obtained by considering a Ty projection instead. 



2. Heterogeneously-doubly-projected mass bound variables 

Notice that in defining the mass bound variables in 
Table |IV[ we have chosen the second level of projection 
(along Tl) to be performed with the same type of trans- 
verse projection ("T", "V" or "o") which was used to 
project into the transverse plane. Of course, this does 
not have to be the case — and by allowing for differ- 
ent types of transverse projections for T and for Tl, 
one would obtain 18 additional variables with "mixed" 
transverse projections. These heterogeneously-doubly- 
projected variables are listed in Table [V] where the addi- 
tional subindex on _L specifies the type of T± projection 



as being of the "T" 



"V" or 



type. As usual, the 



Early partition 


Hedged partition 


Late partition 


MjVT±v , MjVT±<, 
Misio±j- , MjVo±v 


Mt]v±v, Mxivxo 


Mxxviv, Mt±,n 
AfoXyjv, MoXviv 



TABLE V: The 18 additional heterogeneously-doubly- 
projected transverse mass variables for each A*', where the 
additional subindex on _L specifies the type of T± projection 
as being of the "T" , "V" or "o" type. As was the case in 
Tables |III| and |IV[ "partition" implies the combined opera- 
tion of partitioning the objects and agglomerating them by 
summation into composite objects. 



sequence of indices in both Tables |IV] and |V] represents 
the order in which the operations are to be performed. 

For example, A^foAr^T means 

• project all objects using the massless 'o' projection, 
then 

• partition and agglomerate into N parents, then 

• project into the direction perpendendicular to ut 
using the mass- preserving 'T' projection, then, as 
ever, 

• minimize over all values of the unknown momenta 
that satisfy the constraints. 

Interestingly, most of the "_Lo" heterogeneously- 
doubly-projected variables turn out to be related to each 
other and to the corresponding homogeneously-doubly- 
projected variables from Table [TV] For example: 





— Mno±o 


- Mntu 


— MNo±r 


MoA^N 


= Mo^^N 






MoN± 


= MoN±o 




= MNo±r 



M 



TAf_L, 



= M 



.,(150) 
(151) 
,,(152) 



where the last line ( |152[ ) follows from ( |149[ ) and ( |150[ ). 
The one remaining variable Mon±t rather similar to 



M- 



TAf_LT 



M- 



TN±, 



since the difference between them 



may arise only due to nonzero masses {Mi ^ 0) of the 
individual visible particles. 



Another set of 18 additional variables can be trivially obtained 
from TablelVlby considering a Ty type of projection at the second 
level insteac 



IX. PROPERTIES OF THE MASS-BOUND 
VARIABLES 

We should stress that proliferating the number of kine- 
matic variables in the literature is certainly not among 
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the goals of this paper — on the contrary, we empha- 
sise that these variables are different implementations of 
the general principle described in Section [VII A[ What's 
more we will soon begin to reveal further connections of 
these variables to each other (in Section IX C ) and to ex- 
isting proposals (in Section |x|. But before we proceed, 
perhaps now is a good time to summarize what we have 
accomplished so far. 

The previous discussion has hopefully convinced the 
reader that, once the decision on the targeted event 
topology (Fig. [5]) is made, the choice of relevant invari- 
ant mass variables is straightforward and rather unam- 
biguous. Following the general recipe outlined in Sec- 
tion |VII B| one is able to overcome the two main obstacles 
in any analysis involving missing momentum: 

• The fact that the momenta of the invisible particles 
are unknown. To construct a bound, this prob- 
lem is solved by performing a minimization over all 
possible values of the invisible momenta, consistent 
with the measured ^j,. The minimization fixes the 
values of the invisible momenta (e.g. as in (95) and 



(96l) and from that point on, one works with fully 
specified kinematics in the event. Of course, the 
momenta found in the process of minimization, are 
not equal to the actual momenta of the invisible 
particles in the event, although in some cases they 
can be close, see [TO] . 



• The fact that the total number and the masses of 
the invisible particles are unknown. This problem 
is also resolved through the minimization — as we 
have seen in the explicit = 1 examples discussed 
in Section pyil C[ the mass bound resulting from the 
minimization turned out to be a function which de- 
pends only on a set of N summed mass parameters 
(931 or a set of A'' 3-speeds (118), and is insensitive 



to the number of invisible particles or to the fine 
structure of the individual masses Ma or 3-speeds 
Va connected to parent Pa- We set out a general 
proof for general N in Section p^X A| below. 

It should be recognized that for any practical applica- 
tion, there is no need to consider every one of the vari- 



ables in Tables III and IV since some will be better suited 



than others to the particular task at hand. 

For example, we have seen that the "V" -projected 
quantities assume knowledge of the pr and the speed 
of the particle, but leave the mass and Pz undetermined. 
This means that all of the "V" -projected variables in Ta- 
bles [III| and |IV| should be considered appropriate only 
for experimental situations in which the px and speed 
of the particles are known, but nothing is known about 
their masses or longitudinal momenta. Such situations 
may exist - for example if px can be determined from 
the particles' bending radii in a solenoidal magnetic field 
and speed can be inferred from time-of-flight information 
or from the characteristic angle of any emitted Cerenkov 
radiation. However such cases are the exception, rather 



than the norm in current experiments.^*^ In what follows 
we shall therefore give greater attention to the remaining 
three classes of variables: "unprojected" , "T"-projected, 
and "o" -projected. 

We shall denote a generic mass-bound variable as Mjr, 
where the composite index T is made up from (any num- 
ber of) objects taken from the set {A^, T, V, o}. There are 
seven such possibilities^^: 



J- e {n,nt,tn,No,oN,nw,vn} . 



(153) 



For later convenience, we also introduce the generic no- 
tation J^-y for the T-projected variables: 

J"t e {AfT,TAf} (154) 

J-\y for their two "V" -projected counterparts: 

T\/ e {wv, vjv} , (155) 

and To for the two "o"-projected equivalents: 

e {Jvo, oAt} . (156) 

The large multiplicity is partially due to the different 
possible ways to transversify the energy-momenta of the 
composite daughter particles whose masses and 
are typically nonzero. First, one can choose whether or 
not to project, and then those projections can be of type 
T, V or o (see Section III). In addition, as emphasized 
in Section [Vj the operations of partitioning into compos- 
ite particles and transversifying do not commute, so that 
in general we obtain non-equivalent variables simply by 
switching the order of those operations. As illustrated in 
Tables |III|V] we use the ordering of indices on each vari- 
able (from left to right) to indicate the order of the cor- 
responding operations. For example M2T means that we 
add the 1-1-3 vectors first to form two composite visible 
daughter particles and transversely project later, while 
Mt2 implies the opposite — make a T transverse pro- 
jection before forming the composite daughter particles. 



A. Dependence of mass bounds Afjr on l^I, y, etc. 

We have stated that the dependence of the mass bound 
variables, Mjr, on parameters of the hypothesis is always 
confined to a set of N parameters contained within iVl or 
V" etc. We have not yet proved this statement for general 
values of A^, or indicated whether we can demonstrate 



Even in those cases, it is usual that pT is determined from a mea- 
sured track, so one would also expect to be able to reconstruct 
the polar angle 9, from that track, which would permit and 
hence the full l-|-3-vector of the particle to be determined. There 
would be no need to then restrict oneself to the subset of that 
information held by the corresponding "V" -projected quantities. 
The number of possibilities increases to 17 if one allows a second 
level of projection, as discussed in Section [VIII[ 
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this to be true for other classes of projection not already 
discussed. All we have proved, so far, are the following 
statements, which are specific to = 1 and consider at 
most one projection: 



That Ml depends only on lyl {lyli} (sec ( |97| ) 



• That MiT depends only on 1^1 = {iVli} (see ( |106[ )) 



That Mti depends only on 1^1 = {yii} (see ( |111[ )) 



• That Mvi depends only on ^ = {^1} (see ( [T20| ) 



• That Miv depends only on ^ = {^1} (sec ( [T28| ) 

• That Mio depends on no hypothesis parameters'^, 

• That Moi depends on no hypothesis parameters. 

We now seek to generalize the proofs of the above to 
all other values of N, in a manner that does not make 
specific requirements on T. Specifically, we would like to 
prove that: 



jr depends only on $^ = , . . . , (157) 



where $jr is a set of N parameters, of which there is 
one ($^) for each parent Pa, and where the nature of 

depends on the type of projection in J- (which may 
be arbitrary), and on a, but not on the number N of 
parents in total. In particular, we have already seen to 
expect $jr = {iV^i} when T e {i,it,ti}, and to expect 
$jr = {V^i} when T E {iv.vi}, and we are now seek- 
ing to generalize these to results like "$jr — |]yi]^,]^2} 
when J" e {2,2T,T2}" or "^^ = {^1,^2,1^3} when 
J- e {3v,v3}"etc. 

What we will actually succeed in proving is the 
marginally less general statement that: 



depends only on $^ = , , . . . , 
when N = n, provided that Mj- depends only (158) 
on $jr = I when N = 1 



which reminds us that the generality of the desired (but 
unattainable) result (157) is constrained (for any partic- 
ular projection J-) by the need to prove the result for 
the = 1 case. In other words, though the proof of 
( |158 1 found below will be valid for any projection, the 
desired result (1571 will only be true for projections that 



experience simplification in the A = 1 case. 



Note that depending on "no hypothesis parameters" is a special 
case of depending on a very dull set of parameters = {0i} 
which contain no information. 



The proof of ( 158 ) is astonishingly simple. Consider 



an arbitrary mass-bound variable 



Mjr = min 



max [Ma (5f )] 



(159) 



where each parent mass Ma depends on a corresponding 
set of invisible parameters Sf 



(160) 



(compare this to the analogous relations (481 and (50)). 
Let us now perform the minimization in (159) in two 



steps. In the first step, for each parent, we hold the sum 
of the invisible daughters' momenta constant, and min- 
imize over the internal partition of invisible momentum 
between those daughters: 



Mjr — min 

E qaT='pT 



mm 

ieXa 



max [Ma (5f )] 



(161) 

Since the internal partitions over the invisible momenta 
are done independently for each parent, (161) can be 
equivalently rewritten as 



Mjr — min 

Yl qaT=-pT 



max 



min^ [Ma{S^)] 
2^ qiT=qaT 



(162) 

Now we use the assertion that for A = 1 (i.e. for any indi- 
vidual parent) the minimization over internal partitions 

yields a function of a single parameter $^ as opposed to 
the whole set of parameters S^: 



nin _ [Ma (5f )] ^ Ma (^f) • 

qiT=qaT ^ ' 



mm 

Y. qiT=q 



(163) 



Substituting ( |163 ) into ( 162 ), we obtain the desired result 



= min 

Y qaT=pT 



Ma[$^a 



(164) 



which makes it obvious that Mjr can only be a function 
of the set of parameters $jr = ^$1, . . . , ^^|. 

B. Parental masses and upper kinematic endpoints 



By construction, the mass-bound variables of Table III 
are designed to provide an event-by-event lower bound on 
the true invariant mass M^'^^ of the heaviest parent 



M^^'^ = max [Ma] 



(165) 



Such bounding properties are contingent on us being able 
to make appropriate choices when analyzing the events. 
We initially restrict our discussion of the bound to the 
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case where the set of momentum configurations A per- 
mitted under our assumptions is equal to the set of mo- 
menta £ sampled by nature. We observe that to ensure 
A = £ we must (a) correctly reconstruct the event topol- 
ogy (the number of parents, the number and types of 
daughters, and the association of daughters to parents) 
and (b) employ the true values of the parameters $ used 
in constructing of any nijr variable — i.e. for J-\/, 

and lyr*"^"" for J^y or J^jy. Under A — £ conditions all 
Mjr variables are designed to return values smaller than 
the mass of the heaviest parent 



(166) 



From eq. (166 1, it follows directly that, if we were to 



consider the differential distribution of the same variable 
Mjr over all events, the upper kinematic endpoint M- 
of this distribution also satisfies 



max 
J' 



M^"^"" = max [Mjr] 

all events 



(A = £), (167) 



where we make explicit the requirement that the true 
values of the $ parameters are used. 

There remains the important question as to the cir- 



cumstances under which the inequalities in ( 166 ) and 
( 167 1 are saturated - i.e. the conditions for which a mea- 
surement of the Mjr kinematic endpoint will provide a 
determination of (rather than simply a lower bound on) 
the largest parent mass M-p^^. 

We observe that when A — £: (i) that for any se- 
lected event e (z £ the minimization picks out some non- 
empty subset of momenta if^ that satisfy the global min- 
imum; (ii) that C A] (iii) we may define for conve- 
nience Mj^'^^{c) = maxag-p [Maj^{c)] for any configura- 
tion c € A; (iv) that since we minimize over any unknown 
momentum components, the value of Ai']^^^ evaluated 
for some minimum configuration k € K cannot exceed 
the value that would be obtained elsewhere in A (and 
therefore m £); (v) that projections T e {T,V,o} do 
not increase the invariant mass (33 1; (vi) that Mj- can 



therefore not exceed the largest parent's invariant mass 
since 



M 



M 



{A 



£). 
(168) 

The necessary and sufficient condition for saturation 
of (167) is therefore that there exist some event e for 



which two inequalities in ( 168 1 simultaneously become 
equalities. 

maxag-p [M.a{K^)\ it follows from 



Given that Mf""" 



( 168 1 that a necessary condition for saturation is that 



3(ee5,ae7')[Xa(e) = M^-T 



(169) 



There are cases for which ( 169 1 is not satisfied, and for 
which the inequality in (167) must therefore remain un- 



projection discards all previous infor- 
mation about the mass of the 1-1-3 vector being projected, 
and so calculation of Mqn will return the same value that 



would be obtained if one were to set both Mi = (Vi G V) 
and lyla = (Va € 7^). If all events contain massive invis- 
ibles (or indeed massive visibles) in all daughter sets — 
as would be the case for models like _R-parity conserving 
supersymmetry and UED — then ( 169 1 cannot be true 
for any Aioa and so MoAf can only bound from below, 
rather than determine, the mass of the heaviest parent. 
In Appendix |B 3| we prove the event-by-event inequality 
Mmo < MoN, so any bound that is unsaturated for Mqn 
must also be unsaturated for M^o- 

Compared to the J-q variables, the J-t and J- — N 
variables are subject to less stringent conditions for sat- 
uration, because they retain mass information during the 
process of (absence of) projection. Some of the necessary 
conditions can be inferred from the results of in App. |B 1| 
and |B 2| As an example of these less-stringent conditions, 
if \Ta\ 7^ OVa then a necessary condition for saturation 
will include the existence of events e € £ with vanishing 
relative rapidity between and for some a ^ V. 



As discussed in Section VII C some of the widely used 



collider variables like hx, and p^p belong to the J^o 
class of mass bound variables, and as such can generally 
only place a lower bound on the parent mass if nature 
produces heavy invisibles. In order to really measure 
the mass scale of the new particles when massive invisi- 
ble particles are pervasive, one must work with variables 
which retain the dependence on the missing mass param- 
eters and therefore belong to either the M^r or the Mjr^ 
class of mass-bound variables. Other than the full 1-1-3 
dim. invariant mass, other common examples of such 
variables include V^min [50] (discussed below in Sec. X C ) 
the transverse mass uit in the form [52" that accounts 



for the mass of all daughters (sh own i n Sec. XD I and the 
'stransverse mass' mT2 [Z] (Sec. XH| ). 

Before concluding this subsection, we note that the Aljr 
variables are still useful even when the true values of the 
mass IVl (or speed Y) parameters are not known. The 
most conservative procedure in these situations of uncer- 
tainty is to minimize Mjr over the complete physically 
relevant range of any unknown parameter. This leads to 
^ for Mat, Mjr^, and tof^^l for Mjr^. The 
resulting, conservative, Mjr variables still provide lower 
bounds on the mass of the heaviest parent - though those 
bounds will generally not be saturated. 

A more sophisticated treatment is also possible. For 
example if the "true" value of the summed- mass param- 
eter set for the calculation of MnO^) or Mjr^ (lyl) 
were not known — then one could still view the set of 
endpoint measurements for all possible values of lyl^, as 
one constraint among the A^-f-l unknowns M^^^ and iVl^j, 
(a = 1,2,..., N). Not only is this valuable information 
on its own, the derived functional relationship ^''(lyl) 
is, in addition, often sufficient for determining the indi- 
vidual mass parameters l^Ia. The function M^'^^{l)ll), 
when viewed as an A'^-dimensional hyper-surface in the 
(A^ + l)-dimensional space spanned by M^^^ and lyl, ex- 
hibits certain ridge or crease features, which commonly 
originate from the point marking the set of true values 
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of yi [M1I2S]- (The one-dimensional version of this 
phenomenon was originally discussed in [TU1 - H51 1^ and 
is known as the mT2 "kink". Also see [53] for algebraic 
singularity in relation to the kink.) 



C. Relations among the mass-bound variables 

Some of the variables in Table IIIII are related to each 
other, either in general^'^ 



M. 



AfoJ. 



AfT_Lc 



No^_- 



= M- 



TAf_Lo, 



o_LAf 



oJ_TAf 



= M- 



(170) 
(171) 



(172) 



or under some special circumstances, e.g. massless par- 
ticles: 



MTiv({¥la = 0},{M, -0}) 
Mtw^({I^, = 0},{M, = 0}) 



oN, 



oN±^ 



(173) 
(174) 



Given such exact identities like (170), the reader may 



wonder why we even bothered to introduce separately 
variables like M^j- and M^. In our view, such redun- 
dancy is a virtue, since it offers deeper intuitive under- 
standing of these kinematic variables, and allows one 
to think about the same fundamental quantity in dif- 
ferent contexts, e.g. in (l4'3)-dimensions or in (H-2)- 
dimensions. 

We additionally find (see proof terminating in ( |B53 ) 
in Appendix B 3 ) that the mass-bounds from Table III 



obey a hierarchy: 



Mm = Matt > Mtn > Mon > M, 



I No 



(175) 



Similarly, the doubly-projected mass-bounds from Ta- 
ble [IV] obey the hierarchy 

Matt^L > AfTWi > Mt^LAT > AfoiAT > MoATi = MaTo^L. 

(176) 

From these hierarchies, it becomes apparent that the T- 
projected, late-partitioned variables bear a cost associ- 
ated with the insensitivity to the longitudinal momenta. 
By dropping this information we necessarily weaken the 
bound relative to the early-partitioned versions. Inter- 
estingly enough, the order of projection and partition 
has the opposite effect with the o-projection, since both 
longitudinal and transverse information is contained in 
the masses of the agglomerates, and hence by throwing 
away the masses at a later stage, we in fact throw away 
maximal information and are forced to produce the worst 
possible bound! 



23 Previously in ( |107| we already encountered the N = 1 version of 
eq. | |170| l. The genera l proo f for arbitrary A'^ is provided in the 



appendix in equation ( B49 i 



X. CONNECTIONS TO OTHER VARIABLES IN 
THE LITERATURE 



The existing literature is abundant with a number 
of (transverse) invariant mass variables which were sug- 
gested (at various times and for a variety of reasons) for 
the study of missing momentum event topologies (see 
[T] for a recent review). At the same time, the mass- 
bound variables which we defined earlier in Table |III[ 
were meant to be very general, since they target the 
rather generic event topology of Fig. [5] and are intended 
to have as few hidden assumptions as possible. It follows 
that we should be able to correlate the most useful mass- 
scale variables in th e lite rature to one of our mass-bound 
variables from Table III The purpose of this subsection 



is to demonstrate that this is indeed the case. 



A. Missing transverse momentum 

The defining feature of any "missing particle" event 
is the presence of missing momentum (more precisely, 
missing transverse momentum) This is due to the 
production and escape of a certain number of "invisible" 
particles, which are either sterile, or very weakly inter- 
acting, so that they are not seen in the detector. The 
distribution^^ is perhaps the most widely studied dis- 
tribution in relation to new physics searches, especially 
in models with WIMP dark matter candidates like su- 
persymmetry, UED and so on. Eq. ( 
correlate the variable to our il/io varia 



allows us to 
jle as 



(177) 



We see that as Alio is defined more and more inclu- 
sively, it eventually becomes equal to twice the missing 
transverse momentum. Thus in the case of a singly pro- 
duced parent, eq. (1771 allows us to interpret the usual 



variable (more precisely, the variable 2|)^) as a suit- 
ably constructed (in the Mio sense) transverse invariant 
mass of the parent (see also the discussion at the end of 
Sec. VII C 7). In accordance with (1661, in the — > 
limit the upper kinematic endpoint of the distribu- 
tion gives a lower bound on the parent mass in events 
interpreted as single-parent (TV = 1) production. 



A corollary from this statement is that invariant mass variables 
which make similar sorts of assumptions but do not fit into the 
classification of Table |III| are often both poorly motivated and 
sub-optimal. 

The missing transverse momentum is often labelled called "miss- 
ing transverse energy" and labelled fix or _E™"^'^ in experimental 
papers. As previously discussed we prefer to recognize the im- 
porant distinction between energy and momentum, so use the 
symbol 
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Effective mass rrir. 



The "effective mass" variable defined in (84) can be 



also directly related to one of our variables, namely the 
late- partitio ned, "©"-projected variable Mpi di scussed in 

reads 



Sec. VII C 6 The previously derived eq. ( 135c 



^4'i = mis - Ut- 
Therefore, we obtain the correspondence 

Mol > TOoff, 



(178) 



(179) 



allowing us to interpret rrics as a suitably constructed 
(in the Moi sense) transverse invariant mass of a singly 
produced, semi-invisibly decaying parent (see also the 
discussion at the end of Sec. 



VIIC6I 



The comparison between eq. ( 177 ) and eq. ( 179 1 rather 
nicely illustrates the main point of Sec. [Vj that when it 
comes to transverse projections and forming composite 
particles, performing these operations in different order 
yields different results. In the case at hand, when forming 
composite particles before the "o" transverse projection, 
one obtains 2^^, while by forming composite particles 
after the "o" transverse projection, one obtains rries. 



C. Florida V^min and v'iLin ' variables 



As already seen in eq. (100), in the special case of 



= 1, the unprojected mass-bound variable Mjv = Mi 
is nothing but the subsystem V^l-^u^ ' variable from |46j : 



Mi(lyli) 



-(sub) 



(I^l) 



Ml 



ml 



1/2 



(180) 



Restricting to events with vanishing upstream momen- 



tum (ut 



0) , one gets the inclusive 

V^„ii„(¥[i) 



variable from 



lim Mi(]y[i) 



As advocated in Refs. [HI [5^ , practical applications of 
Mi(]Vli) need not be limited to events in which the ac- 
tual number of parents was A^ = 1. The work of [IHl 150] 
showed that in events with N — 2, the peak in the 
Mi(]yl]^) distribution is correlated with the parent mass 
threshold J2a=i 2 -^^fa ; even if the two parent particles 
Pi and P2 are different. 

Note that the mathematical identity ( 170 1 also allows 
us to write 

v^^,''\l^i) = Mit(^i), (182) 
Vs^inC^i) = lim MitC^i), (183) 



relating the A/imin ^^'^ V^min variables to the transverse 
invariant mass quantity Mi-j, which is simply the total 
transverse invariant mass in the event (after accounting 
for the potential presence of any transverse upstream mo- 
mentum ut). 



D. 



Transverse mass 



Perhaps the most popular variable which specifically 
targets a semi-invisibly decaying resonance, is the trans- 
verse mass rriTev, which, as suggested by our notation, 
was first applied in searches for a leptonically decaying 
M^-boson (see, e.g. |B]): 

^Tev = (eeT + ^utY " {PeT + QutT (184a) 
« 2{\p,T\\quT\~PeT-q.T), (184b) 

where PeT (<7eT) is the transverse momentum of the lep- 
ton (neutrino), and in the second line one makes the ap- 
proximation that the lepton and the neutrino are ap- 
proximately massless. Assuming that the W boson is 
produced singly, with zero recoil (i.e. ut = in our lan- 
guage), the neutrino transverse momentum q^T can be 
identified with the measured missing transverse momen- 
tum and (184b) becomes 



'^Teu ~ '^PeTfj, (1 " COS <j)eu) , (185) 

where ipev is the measured opening angle between the 
transverse vectors PeT and pj, . 

In this simple example of a VF-decay, the two daughter 
particles are massless, but the same idea can be easily 
generalized to the case of massive daughters as [5^ 



{M,,M,) ^ M'^+Mt+2{e,Te,T-PeT-q.T), (186) 

where Mg and are the electron and neutrino masses, 
respectively, and 



eeT 



M3 



-PeT: 



(187a) 



(187b) 



Now let us obtain these results with our formalism. In 
general, we have a singly produced (A^ = 1) parent reso- 
nance, which decays to a single (A^y = 1) visible daughter 
particle and a single (Ax = 1) invisible daughter parti- 
cle. Since there is only one particle in each daughter set, 
Vi — {e} and X\ = {u}, there is no need to form com- 
posite particles, so the order of the operations becomes 
unimportant. However, if the daughter particles are mas- 
sive, the two different types of transverse projections give 
two different versions of the transverse mass variable: 

MiV(Af.) = M^i(M,) 



M2 



■p!t 



M?. 



Ml = M^i = {p,T + q.T) 



(188) 
(189) 
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Here eq. (|188| ) follows simply from the general formulas 
( |106| ) or ( |111| with the identifications Mi = Mi = Me, 
PeT, lyli = ^^1/ an d frp — q^ T- Sim ilarly, 



Pit = Pit 



eq. (189) is obtained from either (135a) or (142a 



Now it is trivial to eliminate ut using the transverse 
momentum relation (51) ut = —PeT ~ (fi^T, and show 



that eq. (188) is equivalent to (186): 



Mit{M,) = Mti(M,) = mTe.(Me,M,) 



while eq. (189) is equivalent to (184b): 



Mio = Moi = niTeAMe = 0, = 0). 



(190) 



(191) 



E. Cluster transverse mass variables 

Next we consider a couple of more complicated single 
resonance processes. The first example is /i — !> ZZ — > 
e~^e~vi> where each Z-boson is assumed to be on-shell, 
one decaying invisibly, the other decaying visibly to a 
pair of leptons. For this particular scenario, Ref. [54] 
suggested the cluster transverse mass variable 



Mt,ZZ = (Et^Zi + Et,Z2) - [PT^Zi +PT,Z2 



[PT.e+e 



(192a) 
(192b) 



Note that the 1+3-dimenstional invariant mass of the 
visible and that of the invisible systems have each been 
constrained to be equal to A/f . These two constraints 
reflect the on-shell hypothesis we have chosen to assume 
for each of the two Z bosons. 

Once again, we can obtain this variable from our Mit 
or Mti. In analogy to the case of rriTei^, we have a single 
parent resonance, the Higgs boson h, decaying to a single 
massive visible daughter, the first Z boson and a single 
massive invisible particle, the other Z-boson. This cor- 
responds to = 1, iVv = 1, Nx = 1, Vi = {Z ^- e+e"} 
and Ii = {Z — > vv}. Correspondingly, we identify 
Ml = Ml = Mz, Pit = Pit - PTe+e-, ¥li - Mz 



and f)rp = qT.vu- Then (106) and (111) simply give 



Ml 



'PT,e+e- 



M| 



(193) 



which is equivalent to (192b) in light of the momentum 
conservation relation ut = —PT,e+e- ~1>t- Thus we have 
proved 



Mit(Mz) = MTiiMz) - Mt^zz- 



(194) 



Another interesting example is provided by the process 
h — W'^W~ — e^e~i/P, for which Ref. [SS] proposed the 
cluster transverse mass variable 



M, 



c,ww 



= JM: 



'Pt,, 



-{PT,e+e- +#t) 
(195) 

Here the two leptons are clustered together (even though 
they originate from different M^-bosons, they have a com- 
mon parent in h) and their total transverse momentum is 
PT.e+e- ■ The definition ( 195 ) is similar to ( 192b ), the dif- 



ference now being that the two leptons are not correlated, 
and their invariant mass does not have to be consistent 
with Mz- In addition, the invisible mass parameter l^li 
is now set to zero (as opposed to Mz), because the in- 
visible particles (the two neutrinos) are massless. 

The cluster variable ( 195 ) can be readily obtained from 
Mnj- with the following interpretation: = 1, N\> = 2, 
Nx = 2, Vi = {e~^,e~} and Ii = {v^D}. Correspond- 
ingly we identify Mi = M^+e-, Pit = PT.e +e- and 
]y[i ^ 2M„ = 0. Then the general formula ( |l06[ ) reduces 
to 



'PT,e+e- 



M^(0) = 



which is the same as ( |195| ), so that 

Mit(0) = Mc.ww- 



Ux: 



(196) 



(197) 



Notice that in this example we are clustering two visible 
particles, and the order of operations becomes important. 
Therefore, here Mit and M ti in gener al lea d to distinct 
variables, unlike the case of (190) and (194). 



F. The my""^ transverse mass variable 



Concerning the same h — W^W~ — )• e^e^vD exam- 
ple, Ref. [55] advertized the variable (assuming massless 
neutrinos) 



-Pt ■ PT; 



(198) 



which can be rewritten as 



= A^e'+e- +PT.e+e- +fT' PT.e+e- ^ 



We note that if one wishes to relax the assumption of an on- 
shell Z leading to the visible e+e~ system one may do so by 
treating the electron and positron vectors as separate inputs to 
the visible system V. Similarly one may relax the assumption 
that the invisible system is the result of the decay of an on-shell 
Z by treating the neutrinos as independent invisible inputs. 



^'^C'.WW ■ 



T V ^^e+e- + PT.e+e- " ^^T ' PT.e+e- 

s 2 2 

- {PT,e+e- +0t) 

(199) 
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From ( 197 1 and ( 199 1 it now follows that 



(200) 



This connection in fact was the primary motivation for 
introducing the my"° variable in the first place |56) . 



reader is prepared to handle any assumed event topol- 
ogy, and will be able to construct the proper transverse 
invariant mass variable for the case at hand. 



H. Cambridge mT2 variable 



G. 



The m^"|? 



transverse mass variable 



Our final single resonance example will be taken from 
a new physics scenario, namely a generic model with a 
new Z' gauge boson which decays to a SM Higgs boson 
h and a SM Z-boson as Z' — )• hZ — >• bbiyD. For this 
particular topology, Ref. |57j considered the transverse 
mass variable 



Ml 



PTh 



Ml 



(201) 



where Mh (pTh) is the measured invariant mass (trans- 
verse momentum) of the bb jet pair resulting from the 
decay h ^ bb. 

In our language, the event topology Z' hZ bbvv 
corresponds to a single parent, the Z' boson, thus N = 
1. There is a single {N\; = 1) visible daughter particle, 
which is the reconstructed Higgs boson: Vi = {/i — )■ bb). 
There is also a single {Nx = 1) invisible daughter particle, 
which is the invisibly decaying Z-boson: Xi — {Z — >■ vD}. 
Thus we identify Mi ~ Mi = Mh, Pit = Pit — PTh and 

ly^i = Mz. 

Again, we have chosen to make assumptions about the 
l-|-3-dimensional invariant masses of the visible and the 
invisible systems, requiring the former to be equal to Mh, 
and the latter to be equal to Mz, reflecting our assump- 
tions about the decay topology. As before it would be 
possible to independently relax either of both those as- 
sumptions by treating the bb and vD as independent in- 
puts to the visible and invisible systems respectively. 

If we retain the mass-shell constraint for both the h 



boson and the Z boson then ( 106 1 and ( 111 I give 



Ml^{Mz)^M\i{Mz) 



Ml 



PTh 



Ml 



Comparing to (201 ), we see that 



lim Mi-y{Mz) 



lim Mti{Mz) 



m 



TZ' 



(202) 



(203) 



Since 



'TZ' 



was properly defined as a transverse mass 



variable, it is not surprising that it can be obtained as 
a special case of the mass-bounding variables shown in 
Table III The importance of eq. ( 202 ) is that it shows the 



proper way to generalize m^fj? to the case where the Z' 
is produced inclusively, with some non-vanishing UVM 
ut in the event. 

This concludes our discussion of singly produced reso- 
nances. We are hopeful that after all these examples, the 



The variables considered in our previous examples re- 



ferred either to the event as a whole (as in Sees. XA 



X C ) or to the production of a single resonance (as in 
Sees. X D X G I. We now move on to discussing variables 



intended to handle the production of more than one par- 
ent resonance {N > 1). Such cases are very common 
in new physics scenarios, especially if the new model 
contains a dark matter candidate, whose lifetime is pro- 
tected by some discrete symmetry (typically a Z2). In 
such models (e.g. supersymmetry, extra dimensions, lit- 
tle Higgs theories etc.) the main production mechanisms 
usually involve the pair-production of new particles, thus 
the case of iV = 2 has received the most attention so far 
in the literature, although iV = 3, 4, . . . cannot be ruled 
out, and in principle deserve attention as well. 

A popular variable of this type is the Cambridge mT2 
variable defined as [7] 



mT2 



mm 



[max[XiT,7W2T] 



(204) 



where A^it and M.2T are the transverse masses of the 
two parent particles, and the minimization is done over 
all possible partitions of the transverse momenta of the 
invisible particles, consistent with the measured p^,- 

We note that if there is only one visible particle belong- 
ing to each parent, we can immediately identify the Cam- 
bridge variable r7iT2 with both M2T and Mt2 (and even 
with M2 using (108)) since partitioning and projection 
commute for single particles. If, however, we intend to 
apply TOT2 to events in which one or either parent has two 
or more physical daughters (e.g. when doing top quark 
mass measurements in the di-leptonic ti — > bbl^l^vD 
events) then M2T will become inequivalent to Mt2 and 
we should decide which of these is the right thing to use. 
The answer to this question is subtle. The original mT2 
paper T] does not explicitly state how parent momenta 
should be constructed in the event that they have come 
from compound objects, so it is left up to users to decide 
which inputs to supply. It was certainly in the minds 
of the authors of [7] that users ought always to supply 
the maximal amount of trustworthy information to any 
analysis of any kind. In the context of mT2 this maxim 
would imply projecting only after combining the primary 
(1-1-3) momenta of any constituents of parents, provided 
that those constituents could be "trusted" .^^ Only when 



This might include the case where an experiment that records di- 
leptonic top-pairs with good (signed) fe-tagging could allow the 
Z+ to be associated unambiguously with the b and the l~ to be 
associated with the b. 



32 



defined in this manner (i.e. as M2t) can the maximum 
amount of information be squeezed from the variable in 
"clean" events. However, there can be benefits from us- 
ing A-I-J2 (see for example Ref ^) in high- multiplicity or 
inclusive situations in which the indiviual momenta mak- 
ing up each parent have dubiuos provenance or poorly 
measured longitudinal momenta. In such cases, one can 
benefit from using Mt2, even though its end point is less 
sharp for the signal, simply because it is less sensitive to 
longitudinal momenta and momenta at high rapidities. 

One might ask which of the two 771^2 choices — M2T 
or M-J2 — is "better" . Unfortunately, this question does 
not allow a "one size fits all" answer. Each of the two 771^2 
implementations has its unique advantages and disadvan- 
tages. The longitudinal correlations among the visible 
particles which are preserved by M2T result in steeper, 
better defined endpoint structures — see Figs. 11 [a.) and 



12 'a.) below. On the other hand, Mt2 dampens the ef- 
fects of any longitudinal momenta, which would be ben- 
eficial in circumstances where forward jet activity due to 
ISR may be a problem. 

We hope that the current paper will serve as a reminder 
that studies using 171x2 in which either parent is built 
from two or more reconstructed momenta should think 
carefully about the advantages and disadvantages of both 
approaches before choosing the option that is best for 
them. Both versions of 7tt.t2, namely AI2T and Mt2, 
may prove to be useful, and it can be important to make 
the distinction between them. 

In conclusion of this subsection, we highlight the anal- 
ogy between mT2 and V^min: two variables which are 
more closely related than one might think. We have 
shown that the (1-1-3) dimensional version of mT2, to- 
gether with the mathematical identity ( |170| ) implies 

m^^+''\yi) = M2Tm^M2m- (205) 

The second equality here emphasizes that, in spite of the 
transverse index "T", the 171x2 variable is a bona fide 
1-1-3 dimensional quantity. In other words, the appar- 
ent transverse projection in the definition (204) does not 
lead to any loss of useful information.^^ Of course, the 
same cannot be said about the (l-|-2)-dimensional version 



,(1+2) 
'•T2 



M- 



T2- 



Now compare ( 205 1 to the analogous equation follow- 
ing from ( pOl and (|182[) 



-(sub) 



4i,/(I^i)-MiT(^i) = Mi(Iy^i). 



(206) 



We are reminded that 



-(sub) 



and ; 



Jl+3) 



have essentially 



on a mass in (1+3) dimensions as a function of the 
corresponding invisible mass parameters. In the case of 

'' that mass is the center of mass energy of the 
collision since it views the whole collision as a "single 
parent" , while Tnlj,^^^ hypothesizes that the collision was 
a 2 — >■ 2 process, and therefore bounds the mass of the 
heavier of the two outgoing particles. 



I. The doubly projected variables mT2± and 7727-2 | 

The doubly projected variables 777t2_l and 777^211 intro- 
duced in j25j are nothing but the one-dimensional ana- 
logues of the Cambridge variable 7777-2 ( 204 ) , where one 



performs an additional projection on the directions Tl 
and Ty , correspondingly. As was the case for mT2 dis- 
cussed above, the transverse projection should be inter- 
preted in the T — T sense, and there are three possible 
versions of each variable, depending on the order of op- 
erations: 

{Af2T_L(¥l) in 1 + 3 dims; 
AfT2_L(¥[) in 1 + 2 dims; (207) 
AfT_L2(]yi) in 1 + 1 dims, 

and similarly for 7777-211 (1^)- The example considered 
in was inclusive chargino production, where each 
chargino parent decays to a visible lepton and an invisi- 
bly decaying sneutrino. In this case, each visible daugh- 
ter partition Va has a single massless visible particle (a 
lepton) and the distinction between the early partitioned 

version 777^^21^ ■* = M2t_l, the in-between partitioned ver- 



sion 777- 
(1+1) 

''T21. 



(1+2) ^ 
T2_L — 



Mt2_l and the late partitioned version 

™T2J_' = -^TJ_2 does not become manifest. However, in 
more complicated scenarios with multiple visible daugh- 
ter particles, one would in principle obtain different re- 



7^7. 



(1+2) 
T2J_ 



and 777l^2±^ which is why one 



suits from 7n^2±\ 
is advised to carefully define which particular version of 
wt2_l (and similarly for 7777-211) is being used. 



J. Additionally constrained variables 

So far we have been discussing very general variables, 
which target the most general event topology of Fig. [5] 
Notice that we have made very few assumptions on how 
the decays (52) actually take place, and where such as- 



the same physical meaning: they both give a lower bound sumptions have been made (such as in Section IXEI), it is 



This fact is found to be suprising to people who view mT2 and 
similar variables as acting on "projected" quantities. On the 
other hand it is no surprise to those who have always viewed 7777-2 
as a variable insensitive to relative rapidity differences beteen the 
(total) invisible and (total) visible decay products of each parent 
— the line taken in [9]. 



always possible to relax those constraints if desired. Also 
we did not use any additional information which may be 
available from the preliminary studies of other variables 
related to our events, for example the invariant mass dis- 
tributions of the visible daughter collections Va- Armed 
with such additional information, one may in principle 
further constrain the minimization over the unknown mo- 
menta, and obtain new, more specialized versions of our 
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variables. However, the downside is that such additional 
information typically comes at a cost: the need to make 
additional assumptions about the event topology. 

As an example, consider the Oxford M2C variable p^, 
I17j . which is a variant of ^1^2, subject to the following 
additional assumptions: 

1. The two parents Pi and P2 are identical, with a 
priori unknown mass -Mp^ = Mp^ , therefore the 
minimization over invisible momenta is performed 
subject to the additional constraint 



(Pi + Qif = (P2 + Q2 



(208) 



with Pq and given by (531 and (57) correspond- 
ingly. 

2. There is only one and only one invisible particle in 
each invisible daughter set, i.e. that 

|Xi| = 1x21 = 1. 

3. There is more than one visible particle in each vis- 
ible daughter set, i.e. that 

|Vi| = IV2I >2. 

4. There are no intermediate on-shcU resonances, so 
that the decay (52) is effectively (|Va| + |Ia|)-body 
for each a = 1,2. 

Given a large sample of events S that satisfy these con- 
ditions one can study the distribution of the invariant 
mass Mv_^ of the visible particles in each set Va- If 
presented with some sufficiently large sample of events 
S — {e \ e £ A} one can measure the upper bound for the 
Afva distribution which will be found at the mass differ- 
ence between the parent and the single invisible daughter: 



A™'''' = max [AfvJ = Mp„ - lyl^ 



(209) 



One can then reuse this measurement for any event e £ S 
by asserting the constraint 



A" 



(210) 



during the process of minimization over Q^. 

The advantage of such additionally constrained vari- 
ables is that they are clearly better adapted for the 
study of the corresponding class of more restricted event 
topologies. And additional constraints can bring qualita- 
tively new features, including otherwise unobtained upper 
bounds on parental masses [T71 [22] ■ Their disadvantage 
is that they are better adapted (perhaps only suitable) 
for the study of those restricted topologies, and it is not 
clear how to interpret them once some of the assumptions 
hardwired in their definitions cease to be valid. 



K. Other variables 

In the literature one may sometimes encounter vari- 
ables which have the appearance of a (transverse) invari- 
ant mass, but cannot be related to any of our variables in 



Table III As an illustrative example, consider the trans- 



verse mass variable 



#2 



M2 



(211) 



proposed in Ref. [S^ in relation to th e h — ^ W^W" — >■ 
e^e^iyv process discussed in Sec. XE Comparing to the 
definition (1951 of Mc,ww and to the identity (196), we 



see that Mx^y^y can be formally obtained from Mq ww 
Mit(0) with the rather ad hoc replacement 



T^i =0 



^1 = Af, 



(212) 



However, there is no good physics justification for this 
conjecture and as a result, Mx^y^y cannot be related to 



any of the variables in Table III Not surprisingly, subse- 
quent studies [56] found that Mc,ww = Mi-f{0) outper- 
forms Mt„„ ■ 



XI. SIMULATION: PHYSICS EXAMPLES 

In this section we provide an illustration of our pre- 
vious discussion with two specific physics examples from 
the Standard Model: 

• A case with N = 1. Here we consider the inclusive 
(single) production of a SM Higgs boson (mostly 
from giuon fusion), followed by the decay of the 
Higgs to a leptonic T4^-pair: 



pp^h + X 

^ w+w- 



X 



(213) 



where X plays the role of UVM and stands for jets 
from initial state radiation, unclustered hadronic 
energy, etc. In terms of our previous notation, this 
case involves one parent {N = 1), two visible parti- 
cles (iVy ~ 2) and two invisible particles {Nx ~ 2). 

A case with N — 2. Here we consider dilep- 
ton events from inclusive tt pair production, where 
both W's decay leptonically: 



pp^tt + X 

bbw+w- 
bU+r + 



X 



(214) 



X. 



This case corresponds to two parents {N — 2), four 
visible particles (iVy = 4) and two invisible parti- 
cles {Ni = 2). 
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In both of those two cases, the events very closely resem- 
ble the typical SUSY-like events, in which there are two 
missing dark matter particles. Parton-level event sim- 
ulation is performed with PYTHIA at an LHC of 
7 TeV, including the effects from the underlying event 
(using PYTHIA's default model for it). 



A. An N = 1 example: a Higgs resonance 



of ( 168 ) is (see also Fig. nol 



We start with the Higgs production process (2131 for 
a Higgs boson mass Mh = 200 GeV. In the language of 
Fig. [5] the Higgs resonance is treated as the only heavy 
parent particle (N ~ 1) and the event is partitioned as 

This partitioning is pictorially represented in Fig. |8] We 
now concentrate on the five unprojected or singly pro- 
jected variables which are of interest to us, namely Mjr 
with T S {1, IT, Tl, lo, ol}. Their distributions are 
shown in Fig. 11 ^a), and for proper comparison, we use 
the correct value of the missing mass parameter = 
where necessary. In that case, according to the gen- 



eral property (167), all Aljr variables are bounded from 
above by the parent mass, in this case Mf^. For refer- 
Fig, ~ 



ence. 
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a) also shows the Breit-Wigner distribu- 
tion of the Higgs resonance (yellow-shaded histogram). 
Fig. [Tl|a) c onfirms that the distributions obey the bound 
of eg. ( jl67p . Furthermore, it also shows that each of the 
five distributions appears to be saturated - i.e. that each 
has a kinematic endpoint at the value of the Higgs boson 
mass Mh (only a very tiny fraction of events is observed 
to exceed the bound, but this is due to the finite width 
of the Higgs parent). 

We can confirm the endpoint is saturated for each of 
the variables Mj- G {l,lT,Tl,lo,ol} by explicitly con- 
structing an extremal event. We do so under the approx- 
imation that — Ml, = 0. Following the arguments 
of Sec. IIXBI we should construct an extremal event e 
from the subset £ sampled by nature of the total set 
of momentum configurations A that satisfy our general 
= 1, — Mh topology. In the case of the decay 
of interest ( |213[ ) nature obliges us to impose an addi- 
tional on-mass-shell condition for the intermediate 
bosons £ = {aeA\{P^ + Qif=M^ (z G {1, 2})}.29 
An example of an extremal event e € £ that satisfies the 
constraints and that also saturates the two inequalities 



^® We assume the W to have narrow widths. In fact since we can 
construct an extremal event while imposing a strict on-shell re- 
quirement for the intermediate then we can certainly also 
do so when that requirement is relaxed by allowing the W to 
sample from its natural width distribution. 



Pe+ = (El, El, 0,0) 

Pe- - (£^2,^2, 0,0) 

= (£;2,-£^2,0,0) 

Q„ = {Ei,-Ei,0,0), 



where Ei^2 ^^Mfji-M, 



W-, Q 



-9- 



-e- 



FIG. 10: Illustration of an example extremal configuration 
for the N = 1 variables when applied to the h example. 



We note that in the cases of Mio and Moi, the kine- 
matic endpoint coincides with the mass of the parent only 
because the final state objects in this example happened 
to be massless. In more general scenarios with massive 
particles the endpoints of Alio and Moi will provide only 
an unsaturated lower bound on the parent mass, in line 



with ( 33 ) 



Fig. 11 'a) also allows us to compare the different Af 



J" 

distributions to each other. As expected from the general 
property ( |170[ ), the distributions of Mi and M ij (g iven 
in blue) are identical. As discussed in Section |X C and 
shown in eqs. ( 180 ) and ( 182 1, they also coincide with the 

distribution of t he '\/ lLin variable from |46j . Similarly, 
in line with eq. ( |l73[ ), the distributions of M^iC^i = 0) 
and Moi (shown in red) are practically indistinguishable, 
since the lepton masses are so tiny. Notice that this is 
only true when Myi is calculated with lyli = 0, as was 
done here, otherwise the distributions of M-^i and Afoi 
would generally be different. Finally, the distribution of 
Mio (shown in green) is distinct, as this variable is not 
related to any of the others. 

Upon inspection of the shapes of different distributions 
in Fig. [TTJa), one observes that Mi and Mij- appear to 
peak closest to the parent mass Mh, and consequently, 
have the best defined endpoint structures. On the other 
hand, Mio peaks much farther from Mh, and has a rather 
low event population in the vicinity of its endpoint. Fi- 
nally, the case of Mti(0) ~ Afoi represents an interme- 
diate situation — the peak is found in between the peaks 
of All = MiT and Mio; and the endpoint structure is 
more pronounced than the case of Mio, but not as sharp 
as the case of Mi — M^. This is an inevitable conse- 



quence of the hierarchy (B53) among the mass bounds 



which is present in every event. 

Next, in Fig.[Tl|b) we compare the distributions of the 
standard variables mcff (red solid line) and 2|)^ (green 
solid line) to their mass-bound counterparts Moi (red 
dotted line) and Mio (green dotted line). Recall from 



the discussion in Sec. XB (and in particular eq. (179)) 
that Moi is the analogue of TOoS, and in the limit of no 
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FIG. 11: (a) Unit-normalized distribution of the five TV = 1 mass-bound variables Mj^, J- £ {1, IT, Tl, lo, ol} for the inclusive 
Higgs production process h -s- W+W~ -s- + aX a. 7 TeV LHC, with mn = 200 GeV and lyl = 0. The dotted 

(yellow-shaded) histogram gives the true \/I distribution, which in this case is given by the Breit-Wigner h resonance, (b) Unit 
normahzed distributions of the variables rrics and — (solid fines), contrasted with Moi and Mio (dotted fines). 



upstream momentum the two variables become identical. 
This is confirmed in Fig. |ll[ b), which shows rather simi- 
lar distributions for Afoi and TOoff- However, the analogy 
is not perfect and the TOeff distribution is slightly shifted 
to the right. The only'^'^ reason for this effect is the fact 
that we allow for initial state radiation in our sample, 
so that the Higgs parent is typically produced with some 
recoil and ut ^ 0. This is why the the rries distribution 
does not terminate at M^^ but shows a long tail extend- 
ing to moff > Mh- In contrast, the Moi distribution has 
an exact endpoint at M^. 



A similar analysis holds for the other pair of dis- 
tributions (color coded in green) which are shown in 
Fig. JTlTb). As explained in Sec. XA and seen from 



eq. (|177|), the variable Mio is the analogue of 2f^, since 



the two become identical in the limit of no upstream mo- 
mentum {ut — > 0). However, in the presence of upstream 
momentum, the proper behavior (an endpoint located at 
the parent mass) is retained only by the Mio distribution, 
while the 2|ly distribution picks up a long tail extending 
beyond the true value of Mh- 

This concludes our discussion of the = 1 unpro- 
jected and singly projected variables in relation to Higgs 
production (213). We note that one could also apply 
N = 2 variables to this example, this time considering 
the two W bosons as the two heavy parent particles, and 



We checked that when one restricts the plot only to events with 
Ut = 0, the distributions of rriaB and Afoi become identical, as 



partitioning as 

Vi = {£+}, 

The upper kinematic endpoints of the resulting distribu- 
tions will be found at the corresponding parent mass, in 
this case the mass Mw of the VF-bosons. 



B. An N = 2 example: top quark pair production 

As our next example, we consider dilepton events from 
the top quark pair production process (214 1. We assume 



required by eq. ( 179 I 



that the two 6-jets from the top quark decays have been 
tagged, which distinguishes them from QCD jets from 
initial state radiation. Correspondingly, the tagged 6-jets 
will be included among the set of visible particles, while 
any remaining QCD jets will contribute to the UVM cat- 
egory. 

We first reconsider the = 1 variables already studied 
in Sec. |XI A} and show that they can be useful even when 
there are multiple parents in the event. For the purpose 
of constructing A^ = 1 variables, the event is partitioned 
simply as 

Vi - {b,b,i+,r}, 

Fig. [T2|a) displays the distributions of the resulting 
A^ = 1 variables. Those distributions should be con- 
trasted with the true \/ s distribution of the ti pair, which 
is shown in the figure with the yellow-shaded histogram. 
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FIG. 12: The same as Fig. |11[ but for the tt example. The yellow shaded distribution now gives the true invariant mass of 
the tt pair. 



Just Uke in Fig. 11 ^a), we find o nly t hree distinct dis- 
tributions, since Mi = Mn from ( |170| ) and MtiC^i = 
0) ~ Moi from (173 1. The hierarchical ordering of the 
three distributions is the same as in Fig. 11a), the dis- 
tribution of Ml = A/iT being the hardest, and the distri- 
bution of Mio being the softest. Since all of our TV = 1 
variables are defined through minimization, each variable 
provides a lower bound on the true center-of-mass energy 
a/s in the event. As one might have expected, it is the Mi 
(or equivalently, the A/it) variable which offers the best 
(in the sense of being most stringent and meaningful) 

bound. Since Mi and A/it are identical to the VS^^m 
variable, their distribution exhibits the interesting prop- 
erty first noted in |50) in relation to V^niin: that the peak 
of the distribution is located very near the mass threshold 
for producing the two heavy parents, in this case the two 
top quarks. Indeed, notice how the peak in the (blue) 

Ml = AfiT = V^iiin ^ histogram coincides with the onset 
of the (yellow-shaded) true VS distribution. When ap- 
plied to searches for new physics, one can then use the 

peak in the A/i — Mij- — Vsl^^^^ "* distribution as a rough 
estimate of the new physics mass scale [JSl [SD] ■ 

In analogy to Fig. [Tl|b) here we can also perform a 
comparison of the usual variables TOcff and 2p to their 
mass-bound analogues Afoi and Mio- In Fig. |12[b) we 
compare rrics to Moi (in red) and 2^^ to Alio (in green). 
This time the differences are much less pronounced that 
the single resonance case shown in Fig. [lljb). This sug- 
gests that for N = 2 processes, the variable rrioff (2|^y) is 
on an equal footing with A/oi (Afio). 

We remind the reader that the = 1 variables shown 



Thus let us now discuss the N = 2 variables, which have 



in Fig. 12 'a) do not exhibit any upper kinematic end- 
points, since they are being applied to A = 2 events, i.e. 
they have the "wrong" value of A^ and so the bounding 
relations (166) do not apply to any individual parent. 



the correct value of A^ and for which ( 166 ) holds. In the 
case of A = 2, a tt dilepton event is partitioned as 

Vi = {b,£+}, 

V2 = {b,r}, 

^1 = We}, 

This partitioning can be pictorially visualized in Fig. [7] 
Since we are primarily interested in the kinematical ef- 
fects, for this illustrative example we make the simplify- 
ing assumption (unlikely to be realized in any real exper- 
iment) that each lepton can be associated with its sibling 
&-jet. 

The distributions of the corresponding five N — 2 vari- 



ables are shown in Fig. 13 where for illustrative purposes 
we use Monte Carlo truth information to properly as- 



sign the correct 6-jet to each lepton. According to ( 166 1, 
these distributions are bounded from above by the in- 
dividual parent mass, which in this case is the mass of 
the top quark. Correspondingly, in Fig. [13] the reference 
yellow-shaded distribution now shows the (average) top 
quark mass in the event, which follows the familiar Breit- 
Wigner shape (compare to the Higgs resonance shape in 
Fig. [11]). 

As before, we observe three distinct distributions, 
M2 = M2T (in blue), A/T2({l>^a = O}) ~ A/02 (in red) 
and A/20 (in green). All of them exhibit an upper end- 
point less than or equal to the top quark mass A/j, in 
accordance with ((166), but the three shapes are consid- 
erably different. As before, and in agreement with the 
general hierarchy proven in the arguments leading up to 
(B53) for any event, the early partitioned versions A/2 



and A/2T have the steepest endpoint, with the largest 
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FIG. 13: The same as Fig. 12 but for N = 2 variables: (a) the unprojected M2 and the singly projected variables M2T, 
Mt2, Mo2 and M20; and (b) the doubly projected variables M2t± (black), Mt2± (cyan), Mr ±2 (magenta), M2o± (red), Mo2± 
(green) and Mo±2 (blue). The yellow shaded distribution now gives the average top quark mass in the event. In panel (b), 
"T"-projected quantities are denoted with solid lines and are evaluated with l^Ii — ^12 — 0, while "o"-projected quantities are 
denoted with dotted lines. 



fraction of events near the endpoint. 

Again, we can show that we expect the M2, M2T, M-f2, 
Mo2 and M20 bounds to be saturated by explicity con- 
structing an extremal event e G £ that satisfies the on- 
shell constraints of the t and t quarks and bosons. 
An example of such a configuration (see also Fig. 14 1 is 



Pb = Pl = {Eb, pb, 0, 0) 
= Pe- = {Ee, pe, 0, 0) 
Q. = Q^^iE,,~E„0, 0), 



where. 



Ei = 



^b 

Pb = A(A/t, Mb, Mw) 
E, 



Ej = pj 



Mi 



= p e 



a = sinh 



Pb 



Pi = Ali sinh(p — a) 
p* = \(Mw, Mi, 0) 

Mf, 



sinh 



and the two-body momentum function is given by 

V(a2-(& + c)2)(a2-(6-c)2) 



A(a, 6, c) 



2a 



Recall from eq. (205) that the early partitioned 
variables M2 and M2T are equivalent to the (1+3)- 

dimensional version 171^^2 of the Cambridge 771^2 vari- 
able. It appears therefore that for studies like the one 
presented here, where the UVM contributions can be 
safely identified and accounted for, mi^2^'^ is preferable 

(1+2) 

over rn^2 • 

We use the tt example to also illustrate the doubly 
projected variables from Table IV Fig. [iSjjb) shows the 
doubly projected "T" -projections (solid lines) and the 
"o"-projections (dotted lines). All "T"-projected quan- 
tities are evaluated with = Iyl2 = 0. Each type of 
projection can be done in three different ways: early par- 
titioning, M2T± (black) and M2o± (red); late partition- 
ing, Mt±2 (magenta) and Mo±2 (blue); or in-between 
partitioning, Mt2± (cyan) and Mo2± (green). 



Similarly to the result from Fig. [T3[a), Fig. 13 'b) also 
reveals that the early partitioned, "T"-projected variable 
M2T± has the best defined endpoint structure, which 
clearly indicates the value of the parent mass Mj. As 
for the remaining variables, two are identically equal: 



M' 



2o± 



o2±, 



(215) 



t 

_G — 

i/-e--e 

w+ 



t 

■b 

w- 



FIG. 14: Illustration of an example extremal configuration 
for the N = 2 variables when applied to the tt example. 



which is a special case of the general identity ( 149 ), while 
two others are approximately equal: 

MT±2({l^a =0}) «Moi2, (216) 

where a noticeable difference arises only at low values 



due to the finite mass of the 6-quark — see eq. ( 174 1 
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Mass-bound variable 


Existing 
variable 


N = 1 


= 2 




AfTi(]y[i) 


Moi 


Mio 


M2(]Vl) = Af2T(]yi) 


M2T±(]Vl) 


2f^ = 2$T 








ItT — > 










lyil ^ 0,UT ^ 


ut — 










/ 





























/ 


/ 


Me,M^ -i> 


Me,M^ ^ 






Mt,zz{Mz) 


/ 


/ 










Mc,ww 


lyii ^ 














iVii ^ 












m5f|?(Mz) 




MT — S> 










mT2(]yi) 










/ 




mT2x(]yi) 












/ 



TABLE VI: Correspondence between some of the existing variables in the literature, which were discussed in SectionjXj and the 
corresponding mass-bound variables. A checkmark (/) implies an exact equivalence, otherwise the relevant limit ing condition 



is listed. The last variable mT2±{^a) employs the doubly-projected _L construction described in Appendix VIII 



XII. CONCLUSIONS 



The main "result" of this paper is the proposal made in 
Section VII of a general scheme for constructing and cate- 



gorizing the basic invariant mass variables which are best 
suited for the study of missing energy events at hadron 
colliders. As a demonstration of the utility of this gen- 
eral scheme, in Section [X] we showed how a wide vari- 
ety of widely used kinematic variables discussed in the 
literature can be properly accommodated in our frame- 
work. A short summary of this discussion is presented 



in Table VI which exhibits the connections between the 



variables discussed in Section |X] and the corresponding 
mass-bound variables from Tables IIIII and IIVI The table 
reveals that one can give a new meaning to well-known 
variables like and mcs, which were originally intro- 
duced and defined in a way unrelated to any invariant 
mass considerations. Now we see that the same variables 
allow an alternative interpretation in terms of bounds on 
Lorentz invariants of interest as long as one is using the 
"massless" (o) type of projection for the transversifica- 
tion. 

Another lesson from Table [VT| is that depending on the 
specific topology, the same bound may be constructed in 
different ways. A perfect illustration is provided by the 
variable Mi = Afir- As discussed in detail in Sec. |XE[ 
even for the same final state (two leptons and missing 
energy), the variable Mi = Mij- can emerge as differ- 
ing bounds (either M^^zz or Mc^ww) depending on the 
choice of interpretation of the kinematical information. 

But the value of the proposed scheme is not just in 



the accommodation of existing techniques and variables. 
The primary benefit from our approach is that, having 
understood the main principles behind the construction 
of a good invariant mass variable, the reader is now pre- 
pared to tackle almost any event topology, first by realiz- 
ing what are the proper invariant mass variables for the 
case at hand, and second, knowing how to construct and 
calculate those variables. As discussed in Sections [iTlHVH 
there are a number of choices to be made along the way, 
related to the method of transversification, the partition- 
ing of the event, and the exact order in which one takes 
all those operations. The main guiding principle through 
all this is that at the end of the day, one is always going 
to construct a bound on the mass of the heaviest parent. 
In that sense we are extending the principles and meth- 
ods of construction put forth in [7] for 771^2 and [SD] for 

V^min ■ 

As we have seen, many of the generalized mass-bound 
variables are already in use at the LHC and elsewhere, 
but the majority have, for the moment, the status of 
solutions in search of problems. 



Appendix A: Computer libraries offering 
"transverse" energy and mass variables 

Though libraries should be a repository of human 
knowledge, any careful experimentalist will already have 
recognized that the computer libraries which support 
transverse projection methods for Lorentz vectors do not 
always produce the expected behavior. A selection of 
some of the most commonly used libraries and some of 
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Library 


Object 


Method/function name 


ex 


4 


my 


my 


niTi 


ev 


el 


CLHEP 36 


LorentzVector 


mtO 


mt2() 


- 


- 


- 


etO 


et2() 


ROOT g7] 


TLorentzVector 


MtO 


Mt2() 








EtO 


Et2() 


Fast jet [60| 


Pseudojet 


niDerD ( ) 


inperp2 ( ) 








EtO 


Et2() 


PGS [6T] 














v4et(p) 




Oxbridge 
Mt2 [38] 


LorentzVector 

LorentzTransverseVector 

Mt2_332_Calculator 


ETO 
EtO 


ET2 
EtsqO 


LTVO .massO 
mass 


LTVO .masssqO 
masssqO 


— 

mT2_332() 






UCD Mt2 I3H1 


mt2 


Ea, Eb 


Easq, Ebsq 






getjnt2() 






Defining equation in tiiis paper 
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(2051 





TABLE VII: Tlie versions of the transverse variables used in commonly used high-energy physics computer libraries and codes. 
A brief survey of experimental collaborations' software suggests that most follow the conventions of CLHEP. 'LTV is a shorthand 
for the method getLorentzTransverseVector O . 



their methods for calculating transverse variables can be 
found in Table |VII[ In many cases the method of pro- 
jection used (i.e. "T" or "V") is undocumented and 
can only be determined by excavating the implementa- 
tion. What is more, the names of the methods and func- 
tions in some cases produce output very different from 
what the user might expect. The result is that use of 
a plausible-sounding method can land the unwary user 
with a totally unexpected result - for example the CLHEP 
method called mt() returns the T-projected transverse 
energy (ey = \/ M'^ + p^), not the transverse mass they 
might have anticipated. Of course, because of the right- 
hand expression in eq. ([5|, one might fittingly call this 
quantity a "mass" , but in that case the proper nomencla- 
ture should probably be a "longitudinal" mass and not a 
"transverse" mass. 

To the extent that there is agreement on the conven- 
tions, one can see that the most commonly-used libraries 
(root and CLHEP) use the V convention when calculat- 
ing "transverse energy" quantities. The Tevatron and 
LHC experimental collaborations tend to follow the "V" 
conventions when talking about "transverse energy" in 
calorimeters. For analyses where the transverse mass re- 
ally matters, e.g. for W iiy, the (ROOT and CLHEP) 
libraries have no function to return the 'usual' transverse 
mass of Refs PHSl [S3]: rriT must instead be calculated 
explicitly by the user. 



Appendix B: Mass bounds on collections of momenta 

In this section we present derivations of mass bounds 
on collections of arbitrary momenta, which may be repre- 
sented by unprojected vectors and/or vectors transversi- 
fied by any of the projections T, V and o. These cover the 
cases mentioned in |VII[ and justify the representation of 



multibody decays to visible and invisible particles in the 
form of a pair of composite momenta, where all visibles 
are projected identically (if at all) and all invisibles are 
likewise projected identically, though not necessarily by 
the same method as the visibles. 

The question of what goes into the set of momenta from 
which we wish to generate the parental mass bound is not 
a mathematical question at all. However, once that set 
of momenta is formed, the question of how to calculate 
the best bound making maximum use of the information 
contained in that set is entirely mathematical. It is this 
mathematical question that we solve in the this section. 

In essence, we try to answer the following question: 

Given a particular set of vectors, what is the 
greatest possible lower bound that we can 
place on the mass of any parent particle which 
could have have decayed to daughters charac- 
terized by that set? In particular, how does 
that bound depend on the dimensionalities 
and projection-types of the vectors character- 
izing the information about the daughters? 

We shall denote the answer to that question as A^{. . .}, 
where {. . .} is the set of vectors. We do not wish to 
restrict the set to contain only momenta of the same 
type (e.g. only four-momenta). Instead, we permit 
the set, if so desired, to be a heterogeneous mixture 
containing any number of four-momenta, T-momenta, 
V-momenta, o-momenta or 2-momenta. For example, 
M{At',B'',c'f,d'f,efJ^,g^,hT} would denote be the 
greatest possible lower bound on the mass of a particle as- 
sumed to have decayed to (at least) eight daughters, un- 
der the assumption that the only information from which 
we would wish that bound to be constructed were to com- 
prise: the four-momenta of two daughters a and 6; the 
masses and transverse two-momenta of three daughters c. 
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d and e; the three-speed and transverse two-momentum 
of daughter /; and the transverse two- momenta of parti- 
cles g and h.^^ 



1. Parental mass bounds from sets containing any 
two objects 

Before considering parental bounds from arbitrary sets 
of momenta, we shall first consider the bound one obtains 
for each of the ten pair- wise combinations of the various 
types of vectors, i.e.: 



M{A'',b^} MiAf^^br} 

X{a^f,6f} X{a^f,65} M{af,bT} 

M{a'^,b^} M{a^,bT} 

- - - M{aT,bT} 



where the first equality is simply a rephrasing of the 
meaning of as the minimum mass consistent with 

the constraints. The second equality is from the defi- 
nition of the inner product (or physically the definition 
of the mass), and the third equality is from the defi- 
nition of a vector space (physically representing energy- 
momentum conservation). The fourth equality is a state- 
ment that the vectors and are fully specified, 
so the minimization is trivial (no parameters need be 
changed). The hopefully unsurprising outcome, then, is 
that the best lower bound on the parental mass is given 
by the invariant mass of the two daughter momenta A^^ 
and B^'. 



The 



To calculate the bound 



mass bound 



To avoid imposing a physical interpretation on the vec- 
tors (other than that they are momenta), we generally 
work with A's and i?'s, as opposed to the P's and Q's 
used in the main text. The latter carry implications of 
visibility/invisibility that are irrelevant to the considera- 
tions of this section. 

The list above appears to leave out the massless o- 
projection, but this is simply a special case of the T and 
V projections, so the results for o can be derived from the 
other two cases. In its place, we allow for combinations 
of vectors including transverse two-momenta ar^br, in 
which the timelike component is simply unspecified. It 
will be seen that the bounds from combinations involving 
dT,bT simply emerge to be the massless case. 



a. The A4{A'^ , B^} parental mass bound 

We start with a straightforward case, taking care to be 
explicit about the sequence of operations that will also 
be required for the construction of the bound in the less 
trivial cases. The best parental mass bound'^^ given a 
pair of daughter 1-1-3 momenta A^^ and i?^ is given by 



M^{A'',B^'} = min[Af2] 



min[P''P^] 
= min[(A'" + B^)(A^ 
= (A^^ + BniA^ + B^ 
= {A^' + B''f 



(Bl) 



MiAt'^bf} 



(B2) 



we note that 6" contains partial information about 
some l-f3 vector i?'' which projects to fo" under the 
T-projetion, about which the x and y components are 
known, but the z com ponent, b^ is completely unspeci- 
fied. The bound (B2) can therefore be rephrased. 



M^{A^', 6^} = min [(^'^ + B")'^] 



(B3) 



For the minimization we recognise that provided either 
Ma 7^ or la^l 7^ 0, then M is unbounded above as 
bz — > ±00. Provided that we are dealing with particles 
produced with non-zero transverse momentum (which we 
shall assume hereafter), the solution must then be given 
by the local minimum 







_d_ 

db, 

_d_ 
db. 



( 



B^)2 



Ml 



Ml 



2 EaEb -ar-bT ~ a^b^ 



The minimization selects bz/Es 
has equal rapidity to A'^, 

VB = VA- 



a^/EA such that Bf^ 



(B4) 



To calculate the value of the mass bound we recognise 
that, by the definition of the Lorentz transformation, the 
inner product of A^^ and B^ , which we might denote by 
g{A, B) = A^^g^^i^B'^, is invariant under identical Lorentz 
transforms A of both vectors 



Note that it makes no difference wliettier we use /ly instead of 
h" as an input, as tfie information content of each is identical. 
Note that it is simplest to calculate the bound for the squared 
of the parental mass M^{. . .} rather than for the parental mass 
itself M{. . .}. This difference is of no consequence, and so for 
brevity we will talk only of "mass bounds" in the text, ignoring 
the square. 



g{A,B) ^ g{KA,KB). (B5) 

By letting A be a boost along the z-axis corresponding to 
rapidity change — which will then set both rapidities 
to zero, one finds that 

g(A,i?) = eV^^eV''^-aT-6T. (B6) 
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The best lower bound on the parent mass for daughters 
specified by a 1+3 momentum A^^ and a T-projected 1+2 
momentum b" is then given by 



M^{A^',b'^} 



Ml 



Mi 



St ■ b-i 



{UT + bT)"{aT +bT)a 

(af + bfr. 



(B7) 



c. The A4{A'^ jbr} parental mass bound 

The bound on a 1+3 Lorentz vector with a transverse 
two-vector can be found in a similar manner, but the 
1+3 vector which projects to bx is now given by some 
B which has both unknown z component and unknown 
mass. The bound is given by 



M{A,bT}^M{A,B} 



min [A + B]. (B8) 



A similar argument to that which led to ( B4 ) shows that 



the bz component must be such that ys = Ua- The Mb 
minimization selects Mb = 0, so that 



M''{A,bT} = M'j^ + 2{e^r^^\bT\-aT-b' 



(B9) 



so the bound is formed by turning the transverse two- 
momentum into a o-projected 1+2 momentum. Compar- 
ing with ( B7 1, we see that if 5" is made massless AIb = 0, 



then e 



bT I and ( B9 1 is reproduced. 



d. The 7V/f{ay,6y} parental mass bound 

For each of the 1+2 T-projected vectors, the corre- 
sponding set of 1+3 dimensional objects shares the same 
transverse components and inner product (mass) as their 
T-projected counterpart, but has arbitrary z momentum. 
The bound is then given by 

M^{aT,bT} = mmUA + B)^] (BIO) 

a~ ,bz 

This time the minimizations force the rapidities of A and 
B to be equal, but leave the value of that rapidity ~ 
Ub free. Similarly to the previous cases. 



M'^{aT,bT} M[ 



Mi 



CLT ■ bj 



a \2 



(Bll) 



e. The A4{a" jbr} parental mass bound 



The limit on M is given by 

M'^iaT^br}^ min [{A + Bf]. 



(B12) 



The minimizations set the rapidities to be equal yA = Vb 
(but undefined) and Mb — 0. The limit again appears 
in the form, 



M'^iar^br} 



Ml + 2 



{a^+K?. (B13) 



/. The A4.{aT,hT} parental mass bound 

The limit on M for a pair of transverse two-momenta 
is given by 

M^{dTM}= min \{A + Bf]. (B14) 

az,bz,MA,MB 

The z minimizations again set the relative rapidities 
equal but arbitrary yA = ys, and the mass minimiza- 
tions set Ma — Mb ~ 0. 

M'^laTybr} ^ 2 ^IotI |&t| - St ■ 6t) 

^ (a„" + 6^)2. (B15) 



The A4{A'^ ,b"} parental mass bound 



The V projection described in section |IIIB| maps all 
1+3 vectors B^^ with the same transverse momentum 6^ 
and velocity Vb = \b\/EB to the same 1+2 vector by. 
Therefore the longitudinal momentum component bz is 
unspecified and the parental mass bound A^{A'^,6y} is 
given by 

M^{A'',bZ} = min [(A + 5)2] , (B16) 

From equations [22| and [23} we see that we can decompose 
the full (l+3)-dimensional energy and mass 

El = {e^f + {ef)\ (B17) 

Ml = [m^f + (mf )2 . (B18) 

each in terms of a transverse quantity (e^ , m^) and a 
longitudinal quantity (ef = \bz\/VB, mf = \bz\/{ VBlB)) 
with jB denoting the Lorentz factor l/y/l — Vg. 

Using these relations, we can write the Lorentz- 
invariant quantity (A + B)'^ as 



{A^ + Bn'' - Ml + {m^f + bl/{VBjB) 
'Ea 



(B19) 



2 ( b^ + bl - qt ■ br - azb, 



Leaving aside the trivial case of bx = 0, we now at- 
tempt the minimization over bz, requiring 



= A(^p + sM)2 

obz 



(B20) 



= 2 



Ea bz 



Vb ^/yf+bl 
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This gives rise to a quartic in 6^, 

{b^T + bl){h - a f - eHl ^ , (B21) 
where the constants 



Tr2 2 

Ea Vb ll . 



(B22) 
(B23) 



The need to solve this quartic makes the M.'^{A^^ ^ 6"} 
bound intractable in comparison with the similar 
M.^{A'^,h'^} bound. Similar difficulties are encountered 
in the following M.^{ay^ by} case. 

One might guess that the ys = Ua condition result- 
ing from A^^{A^, 6" } bound could represent the correct 
solution, since we are again working with a fully (1+3)- 
dimensional vector combined with a (l-|-2)-dimensional 
vector. But the solution this condition gives for b^ is not 
a root of the quartic in (B21 1. One can show that b'^ = 
implies 



bra 



(B24) 



which when substituted into the LHS of (B21 ) gives 



(a2 - £2^^1)2 



2bTV^JeWi 



which is in general non-zero, i.e. the equal rapidities con- 
dition only minimizes b^ under certain special conditions. 



h. The A^{av,bv} parental mass bound 

In the case of two V-projections, the mass bound is 
given by 



M^{aZ,b^} = min UA + B) 



(B25) 



the minimization of which involves finding and 6^, 
with both Va and Vb being held fixed, such that each of 
az and b^ is a root of a quartic like that in (B21 ). 



In this situation, we can apply the same mass/energy 
decompositions (B17) and (B18l to a^, to get a variant 
of (IbTI), 



+ al ^b\ + 62 



Va 



Vb 



V^Yb 



dx ■ bj 



a,b^ 



Differentiating by az and by bz separately, the mini- 
mization imposes 



= -^M^ + B^)2 
oaz 



(B27) 



2 ( Vbj, + bl az 
[y'a^ + alVAVB 



bz + 



and simultaneously 
d 







db. 



■\/a^ + a2 ( 



(B28) 



+ fo? VaVb 



N ote th at the fraction |a|/|6| appears in both (B27) 
and (B28l, albeit as a reciprocal in the latter. So, we can 



combine the two minimization constraints in the form of 
a quadratic in az and bz'- 



3_ 



This can be solved to give = —cbz, with 



1 or 



(B29) 



(B30) 



Substituting this solution back into ( B26 1 gives a pleas- 
ingly simple result 



imCr + {m^y+bi[2c + 



Va 



Vb 



(B31) 



Since c was chosen to be positive, this expression is 
clearly minimised for 6^ = 0, implying that = as 
well. If wc then make the replacements ot/Va — and 
bx/VB = Cv , we find that the choice az = bz = gives, 
quite simply and in tune with our intuition and inductive 
sense, 



a\2 



i. The A^{ay,bv} parental mass bound 



(B32) 



The bound A4^{a^,by} requires minimization over 
both ttz and bz- 

■^^{ot, = min [(A'' + B^' f] , (B33) 



(B26) with {At" + B^)2 defined as before in pTg]). First the 



minimization over forces the rapidities of A and B to 
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be equal, a^/ Ea 
obtain 



bz/Es- Plugging this into (B19I, we 



+ = Ml + {m^)^ + bl/iVBlB)^ (B34) 

+ 2 (eaEb -dr-bT-bl 



hi — 

\ 'Eb 
= Ml + {m^f + bl/{VBlB? (B35) 

+ 2(^^{El-hl)-aT-bT 

The expression {E'^ — 6^) can be written one of two ways 
- either as (ef )^ or as {e^Y ~ b1/{VBlBY ■ We choose 
the latter, since we have fixed , but if we were to fix 



instead , we would rederive ( B9 ) 



A further simplification is implied by the equal rapidi- 
ties condition, since Ea = e4 Eb/^EI - bl- Using this 
and (B17), we find 



+ = Ml + {m^f + bl/{VBlBf (B36) 



2 e. 



A „B 



bl 



Recognising that 6^.7^ is positive, we see that the bz 
minimization simply gives a-z = bz = 0, and thus 



Q\2 



j. The M{a" ,bT} parental mass hound 



(B37) 



This mass bound is similar to A^{a", 5" } with the ad- 
ditional minimization of Mb ^ 0. 



M'^ia'^M} = min [{A + Bf] 



(B38) 



{m()y + 2{e()bT-aT-bT) (B39) 

-.a I LQ:\2 



(a? + b^y 



(B40) 



2. Arbitrarily large sets of (1+3)-, (1+2)t- and 
2- vectors 



The generalization of (Bl) to an arbitrarily large set 
of fully specified 1-1-3 vectors A = {Af | 1 < i < |^|} is 

M^{A} = min [M^] = min [P^] 

= min[(E,ylf)2] (B4f) 

Note that in the special case of fully specified f-|-3 vec- 
tors, the mass bound for the set is the same as the mass 
bound of the single object formed of the sum of those 
vectors 



m{A} = m{j:,a,}. 



(B42) 



Let us further generalize our results to an arbitrary 
set of (l+3)-vectors A and (1 + 2)T-projected vectors 



Bt = Ifc^T I 1 < j < |St|}. Each of the b'^-y has a (1+3)- 
vector equivalence class B^" for which the z components 
can take any value. Writing B = {Bj \ 1 < j < \Bt\} 
and Bz — {bjz | 1 < j < \l3z\}, we can therefore write the 
mass bound as 



M{A, Bt} = M{A, B} = min 



vector K^+i 



gives a bound M{Kj_^ 



which has free z 



(B43) 

where each of the B'^ has a free z component. The result 
can be found by induction. The bound M{K^, . . . , K^} 
for some set of fully specified (l+3)-vectors is given by 
the sum A^{Ei=i,™ii:f } by ( |B42[ ). Adding a further 1+3 

momentum to that set 
A similar argument 
to that which led to ( B4 1 shows that the rapidity of Kl/^^^ 
must be equal to that of Yii=i^„iKf . With this constraint 
applied i^^+i becomes a fully specified 1+3 vector, so we 
can treat it as one of the known 1+3 vectors and proceed 
with the next 1+2 (T-projected) vector in the set. 

Applying this argument sequentially to the Bj we find 
that 



x{Afo-T}-A^{(s.Af),i?n, 



(B44) 



where each of the Bj has the same rapidity as S^A^. 

Since the set of 1+2 T-projected vectors is isomorphic 
to the set of 1+3 vectors with fixed (but arbitrary) rapid- 
ity under the operations of addition and inner product, 
we can rewrite this bound as 



XM,ST} = ^{(S.Af),(E,6"T)} 



(B45) 



the bound for the summed 1+3 vector (S^Af ) and the 
1+2 T-projected vector (Ej6°y), the explicit formula for 



which is given in (B7). 

We can further extend the argument by allowing 
some other daughters parameterized only by their two- 
momentum to be added to the set, 



M{A,Bt.Ct), 



(B46) 



with Ct = {ckT I 1 < fc < |Ct|}- Each 2- vector Ckr has 
a corresponding equivalence class which can be repre- 
sented by a 1+2 T-projected vector c°-p with unknown 



mass. The arguments which led to (B44) apply equally 



to the CkT, so the corresponding 1+3 vector rapidities 
are set equal to E^A/Xi, but now we have the extra min- 
imization over the masses which fixes = for each 
(or indeed c^-p) . 
Therefore 

X{ASt,Ct} -X{(E,Af),(E,&^T + S,c^J} . (B47) 



Now (B47l has the same form as all the previous bounds. 



but in obtaining the result we have found out something 
non-trivial: one would not get the best bound on M if 
one were simply to replace the set of 2- vectors Ct by their 
sum {Ti}~C}~t)'- One must instead add the corresponding 
massless 1+2 vectors c?„. 
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In principle we could now try to extend our bounds to 
include (arbitrarily large numbers of) V-projected 1+2 
vectors. However we shall not do so for two reasons. 
The first reason is that in collider experiments such as 
the LHC, situations for which V projection is appropri- 
ate are rare. It is only in very unusual cases where we 
might find ourselves knowing just the transverse momen- 
tum components and the size of the three-velocity, but 
not the azimuthal angle 9, the z-momentum, or the mass. 

The second reason we do not pursue the V vectors fur- 
ther is that one ends up with a real mess, as we have 
seen. The most basic pair- wise combination A4{A,b\/} 
requires solution of a quartic equation in ■ Only if one 
is solely interested in combining (T, V, o)-projected vec- 
tors might the expressions be tractable, but the utility of 
such a combination is unclear. 



3. Mass bound hierarchies 

The similarity in the expressions for the mass bounds 
derived in the preceding sections allows for a further ob- 
servation - that as progressively more information is ne- 
glected or unknown, the mass bound is lowered. Intu- 
itively one would expect this, since the absence of hard 
information causes one to have to be progressively more 
conservative, but we can, with little additional work, 
show this explicitly to be the case. 

We set out, t here fore, to prove the hierarchy that was 
seen earlier in (1751. Our proof proceeds in two stages. 



with 



In the first stage we demonstrate the result for the case 
iV = 1, in which the hierarchy becomes: 



Ml = MiT > Mti > Moi > Afio. 



(B48) 



In the second stage we extend this to general N. 

Using the results of the previous section, we can treat 
each of the mass bound variables in terms of the com- 
posite visible and composite invisible objects described 
in section VI B The equality in ( B48 ) then results from 



the definition of Mi as a concrete case of (B7), where A'^ 
represents the visible P^, and 6" the invisible g". Sim- 
ilarly, Mil is just (Bll), where a" a nd b° stand in for 
and g". On comparing (B7) with (Bll), we see that 



they are identical, and hence Mi — A/it. 

For the next statement, Myj > M-^i, we have to con- 
sider the difference between "early" and "late" partition, 
i.e. whether we retain information about the relative lon- 
gitudinal momenta of the visibles. Let our visible com- 
posite of parent Pa be composed of constituents P/', 



ieVa 

(Ea,PaT,Paz) , 



Ea — ^ Ei , 

ieVa 



We form the early-partitioned composite 



PaT 



(eaT,PaT) 



and the late-partitioned composite 



differing only in their energy components 
eaT = V^l - pL 



era 



+ PTa : 



where 



m 



Ta 




Of c ourse, it is established in the preceding sections [B 1 a| 
and Bib that > m^^, since niTa could be con- 
structed by repeated minimizations of over the lon- 
gitudinal momentum components (pz)i- Hence, e^T > 
era- 

If we now define analogous quantities Gar, ^Ta and claT 
for the composite invisible particle, then all the same 
arguments apply. Forming the two mass variables as in 
(fT02| and pXi9|, 



Aff-T- = (eiT + eir)^ - (pit + Qit)^ 

> 

(eTi+^Ti)^ - (PTI +qTi)^ = Af^i . 

Moving next to Afoi, we note that this is simply the 
previous case, with an additional minimization over the 
masses Mi of the constituent particles, which must reduce 
the size of the bound, forcing A/oi < M-fi- 

For the final inequality, we recall the statement due 



to (B47l, that says the bound is weakened (i.e. made 
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smaller) if we base the bound on the sum of the transverse 
two-vectors, rather than promoting them to o-projected 
(l+2)-vectors before summing. The difference is solely 
in the energy component - the late-partitioned p"^ has 
energy component 



eoa = ^ {PT)i , 





= min 




= min 




= min 




= min 




= min 



ieVa. 

whereas the early-partitioned p^^ has energy component 

Goo — PaT ■ 

By the triangle inequality, Goq > Gao, yielding the final 
required result, that 

= (eoi + Goi)^ - (pTi + QTi)^ 

> 

(eio + eio)^ - (Pit + Qit)^ - . 

Armed with this knowledge, we tackle the hierar- 
chy when N > 1. We revisit the definitions of 
Mjv, MnttMtnjMno, and Mon, from section [VII C[ as 

max[A^a(Pa,Qa,Ma)] 

a 

max [M aT (PaT , QaT , P-a) 
a 

max [MTa (PTa , QaTa , P-a 

a 

max [A^ao(Pao, Qao, Aa) 

a 

max [MoaiPoa, Qoa, Ma) 

a 

At first glance, it might seem alarming that we assert 
Ml = MiT, when Mi seems to be built of a fully (1-1-3)- 
dimensional object 7Wi(Pi, Qi, /ii). But in fact, with 
the components {qz)i left free, the minimization will (for 
reasons identical to those in the discussion of early and 
late partitioning) be achieved when all the constituents 
of Qi have equal rapidity to Pi, meaning 

A^l(Pl,Ql,/il) = A^iT(PlT,qiT, Ml)- 

But this should apply to all N, since the only con- 
straint on the invisibles of each parent is on their 
transverse momentum components. That is, for each of 
the N parents, given our inputs we will get 

Ma{Pa,Qa,Pa) = 7W IT (PlT , QlT , Ml ) , 

and therefore we immediately see that 



M^m 



mm 

Yl qiT=pT 

= min 



max[A^Q(P<j,QQ,/ia) 

a 

max {M qT (PaT , QaT , pa 



= Mwt(D^) . (B49) 
Next one might ask whether the successive inequal- 
ities still hold. The very first one follows straightfor- 
wardly. Only in the input vectors to each of the N 



parental mass bounds A^aT do Mjvt and My at differ. 
Furthermore, since the late-partitioned input vectors pTa 
and Qto will have smaller energy components than their 
early-partitioned counterparts PaT and QaT, each of the 
individual parental bounds follows the relation 

Ma\{Pa\-,<\a\,Pa) > A^Ta (PTa , QTa , Aa) , (B50) 

for every possible choice of unprojected inputs P^, Q^. 

To complete the argument, we need to establish that 
the global minimum considering all trial qiT cannot in- 
crease if any or all of the parental bounds decrease. 

The minimization probes the full space of {oiTji sub- 
ject to the constraint that their sum is the missing trans- 
verse momentum vector, with all other parameters hav- 
ing been specified. For the minimization to pick out a 
larger value for A/ttv than for Mjvti we must have 

max[7UTa(PTa,qTa,Ma)] > max [A^^T (PaT , QaT , Ma )] 

a a 

(B51) 

for the same values of {qir} that give the value of Mnt, 



if nowhere else. But we have already established (B50) 



for all a and all inputs. So we are led to the conclusion 

A%T(^)>MTAr(I^) . (B52) 

Actually, we have achieved more than that. The same 
argument holds for the remaining levels of the hierarchy 
involving the o-projection. So we can boldly claim our 
final result and can retire to a well-deserved cuppa 



M 



N 



> 

Mtn 
> 

MoN 
> 

Mno- 



NT 



(B53) 
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